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Introduction 


As you know, one of the main aims of this course is to help you to develop your 
appreciation of the role of mathematics in understanding and predicting the behaviour 
of the real world (as opposed to the world of mathematical theories), and to begin to 
acquire the skill of using mathematics to help to solve real-world problems, which we 
call mathematical modelling. This unit is devoted to that task. 


There was a time when the application of mathematics was restricted more or less to 
physics and engineering; but those days are long gone, and now mathematics turns up 
all over the place. In case you need to be convinced that this is the case, you may like 
to consider the following list of titles of conferences organized in the course of a year by 
the Institute of Mathematics and its Applications. 


Mathematics in Industrial Maintenance 
(topics included component replacement problems, availability and downtime 
modelling, and inventory and manpower problems) 


Business Modelling 
(themes: the mathematical modelling of production and manufacturing, finance, 
human resource management, and marketing) 


The Mathematical Theory of the Dynamics of Biological Systems 
(which covered the epidemiology of infectious diseases, including AIDS, as well as 
the harvesting, control and regulation of biological populations) 


Control Theory 
(economic systems and neural networks were among the topics here) 


The Mathematics of Surfaces 

(which was concerned with the representation of complicated surfaces, and dealt 
with computer-aided design and machining, robot vision, cartography, and even 
applications in television and the film industry) 


Stably Stratified Flows 
(it may not be apparent from the title, but this conference was actually about the 
effects of topography on the weather) 


Mathematics in Food Production, Processing and Preservation 
(covering such topics as freezing and chilling, cooking, and extrusion and injection) 


Aerospace Vehicle Dynamics and Control 
(everything from airships to spacecraft) 


Mathematics in Signal Processing 
(with applications to astronomy, communications, medicine, meteorology, 
radar, etc.) 


Puzzles in Quantum Theory 
(physics has not been completely overwhelmed by newer areas of application) 


All human life is here! Though these diverse uses of mathematics call on different 
aspects of the subject, there is a range of skills which is common to all of them, and 
these are the skills that together constitute the craft of mathematical modelling. 


This unit takes an incremental approach to the teaching of mathematical modelling. In 
the first section you will be introduced to the art of interpreting and analysing a model 
produced by someone else, and in the second you will be asked to carry out this process 
for yourself. In the third section you will practise some of the individual skills of 
modelling. The fourth and fifth sections are devoted to integrating those skills. 
Mathematical modelling can be learnt only by doing it, so there are many exercises in 
this unit: in particular, in Section 5 you are asked to develop a simple mathematical 
model yourself. Sections 1 and 4, however, consist mostly of reading; but I do hope that 
you will read them actively, rather than passively — by which I mean that I hope you 
will try to answer for yourself the many questions that are implicit in the text, and try 
to anticipate the direction that the argument will take. 


Ac=0. 10633 exp ‘i g} 
(1+ L6pq) 


(ph +0.65pr+0.88pm) 
yy-279 


exp {by/Ag+b/(X+2.6) 


a year of glorious follies 
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Two television programmes on modelling will be broadcast during the two weeks that 
are set aside in the course calendar for your work on this unit and the associated 
assignment. Each of these programmes will further illuminate the ideas discussed in the 
unit. Neither requires you to do any specific work either before or after the broadcast, 
and unlike the other programmes for the course, neither introduces 'mathematical 
results which you will need to refer to later. I have therefore not thought it necessary to 
include programme summaries in the unit. 


1 First thoughts about the modelling process 


Some years ago, The Sunday Times magazine for the last Sunday of publication of the 
year appeared with the cover shown here on page 6. Among the other ‘glorious follies’ 
to be found inside were photographs of the Prime Minister with his trousers held up by 
@ paper clip, Prince William falling off a horse, Benetton adverts, and sumo wrestlers. 
On page 3 of the magazine the reader was told that ‘the probability of having a road 
accident (see cover) is only slightly more difficult to calculate than the running total of 
Liz Taylor's husbands’. 


Reading the magazine sections of Sunday newspapers might very well be regarded as 
folly, and not very glorious at that, so perhaps one shouldn't be too surprised when 
mathematical modelling is made into a joke so prominently in one of them. I don’t 
know how you react, but it took me several minutes to recover from my annoyance that 
a so-called quality newspaper should take such a dismissive attitude to what is clearly a 
‘serious investigation of a serious subject. I suppose that I should have learnt over the 
years that mathematics is rarely taken seriously by the media, even though it plays 
such an important role in modern life, 


Once I had got over my irritation, however, | began to think about the formula itself, 
to ponder its significance, and to try to understand what it means. Many of the points 
that I want to raise in this unit came up during my consideration of the formula, so I'd 
like to recount some of my thoughts by way of introduction to the unit. By describing 
my reflections about the formula, I can bring to your attention again some of the points 
about modelling that were first raised (in this course) in Unit 3, and in particular to 
say once more what mathematical modelling involves and why it is important. I should 
make it clear from the outset that the model we are about to consider has been 
produced by professionals at the Transport and Road Research Laboratory, and is more 
complicated and sophisticated than any model you will be asked to develop yourself in 
this course; but I don’t think you will find it difficult to appreciate the main issues it 
raises, 


This section could be described as an exercise in interpreting a mathematical model: 
not an easy task in this case, since there is so little information to go on. There is not a 
great deal for you to do except read; but you will probably find it valuable to think 
about how the points that arise in the discussion apply to other models that you have 
met, earlier in this course or elsewhere. In the next section you will find some extensive 
exercises in which you will be asked to analyse and interpret some mathematical 
models yourself. 


The formula we have to deal with is this: 


A, = 0.006 33 exp{s + g}(1 + 1.6pa)(pp + 0.65p, + 0.88pm, )M°?” exp{bi/Ag + b2/(X + 2.6)}; 


it is a mathematical model of the frequency of traffic accidents. Here already is my first 
point about modelling: a mathematical model is a formula. Perhaps that is a little too 
sweeping, but it does express a basic truth: the aim of mathematical modelling is to 
capture the essentials of some situation or process in the form of one or more equations 
or other mathematical relationships, in order to understand it better, to predict how it 
may develop, or to help one make decisions. 
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The Sunday Times dismissed the formula as a ‘glorious folly’, implying perhaps that it 
has no use. We may reserve judgement on that, but one thing is certain: there is no 
way of telling how useful it may be without being able to say what it means, and The 
Sunday Times has made this task more or less impossible by not revealing what the 
symbols stand for. You and I can see the general structure of the formula: the 
right-hand side consists of a product of several terms, some of which are linear 
functions of certain variables ((1 + 1.6pz), for example), some of which involve 
exponentials, and one of which is a power. The formula ‘enables you to calculate your 
mathematical chance of having a traffic accident next year’. (Why ‘mathematical’ 
chance, I wonder ~ what is an unmathematical chance?) So A, is presumably the 
chance of an individual having a traffic accident during one year. I have deliberately 
changed the wording here: The Sunday Times seems to suggest that the formula is 
personal to you (the reader), and applies only to the year following publication, which 
is hardly likely to be the case. In order for you to calculate your chance of having an 
accident next year, you must substitute those values of the variables s, g, pa and so on 
which are appropriate to you in the year in question. We aren't told what kind of 
accident, incidentally; my guess is that the formula is concerned with motor accidents, 
if only because accidents to pedestrians or cyclists generally seem to be of much less 
interest. to journalists. It also seems likely that the formula deals with drivers rather 
than passengers, and with cars rather than with other vehicles. 


So the formula expresses how an individual's chance of having a car accident depends 
on certain other variables; but what those variables are one cannot tell just by gazing 
at the front cover of The Sunday Times magazine. I became sufficiently interested in 
the question to write to the Transport and Road Research Laboratory for more 
information. But inquisitiveness won't wait for the post (and in any case Christmas 
and the New Year had to pass before I could expect a reply); so I found myself, in the 
meantime, thinking about what factors do affect one’s chances of having a road 
accident. I found it a quite instructive exercise to try to list the things about an 
individual's motoring which would be likely to influence the probability of having an 
accident, and to ask how, qualitatively, they might affect the answer. For example, how 
much one uses a car is obviously important: the more one drives, the more likely one is 
to be involved in an accident. That is worded rather loosely, but the idea could easily 
be made sufficiently precise to make it worthy of representation by a symbol in an 
equation: annual mileage (or, if you prefer, annual kilometrage), You may wish to stop 
reading at the end of this paragraph and spend a few minutes trying to extend the list 
for yourself. It is probably better not to try to guess what the symbols in the formula 
mean, but to think about the problem from scratch. 


I listed the following features, in addition to annual mileage. 
e The driver's age and state of health, 
e The driver's skill, tendency to drive fast, etc. 


© The type of road encountered (it is often claimed that motorways are safer than 
other roads, for example). 


© The age and condition of the car. 
© The driver's gender. 


I don’t claim that this is an exhaustive list, but I do think that all these are relevant 
factors. Whether they could all be defined in such a way that they could appear in a 
formula like the one under consideration is a matter that requires further thought. One 
way of approaching this question is to ask whether and how each factor can be 
quantified. Age clearly can; state of health, and skill, are not so easy to measure 
directly. On average, one’s state of health tends to decline with age (alas), so in the 
interest of simplification it might be as well to concentrate on the single variable age. 
Likewise, to include skill I would have to think of some other, more easily measurable, 
factor to which it is closely related. Gender may appear to cause difficulties: it does 
strain the English language a bit to describe someone’s gender as a variable. 
Remember, however, that the purpose of the exercise is to obtain a general formula 
which can be made to apply to any individual by giving the symbols (or variables) the 
appropriate values. It isn’t too difficult to see how this might work for the person’s 
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gender: one merely has to allow the corresponding symbol to take only two values, say 
0 for a man and 1 for a woman. 


My second point about modelling is this. If the primary aim, in modelling, is to 
describe a situation in terms of an equation or similar mathematical relationship, then 
it is clearly of the greatest importance to decide what variables are going to enter the 
equation. Very often (as here) one will want to investigate one particular variable — in 
this case, the chance of having an accident. It is then sensible to ask oneself what 
factors it might depend on. The answers must, of course, be capable of clear definition 
and measurement, since they are to appear in the equation as variables. One should 
keep at the back of one’s mind the possibility that a variable may have only a limited 
range of allowable values: perhaps as few as two. Incidentally, the variable A, is called 
the dependent variable, because it depends on the others; the other variables are called 
independent variables. 


When I’m thinking about the variables that enter a mathematical model, I make a 
practice of trying to see how, qualitatively, each will affect the outcome. For example, I 
think it is pretty well bound to be the case that older and more skilful drivers are likely 
to be safer. I could express this mathematically by saying that the chance of having an 
accident decreases with increasing age, all other factors being equal; or that the chance 
of having an accident is a decreasing function of age. Another way of representing the 
dependence of the chance of having an accident on age would be to use a graph: I hope 
you can see (in your mind’s eye) what such a graph must look like, in general terms. 


Here is my third point ~ a hint rather than an instruction: you can gain some useful 
initial insight into a model by thinking about the qualitative nature of the relationships 
between the variables. 


So far, I have been describing my speculations while I was waiting to hear from the 
Transport and Road Research Laboratory (TRRL). I did eventually receive a reply to 
my enquiry, and found that I was correct in assuming that the report is concerned with 
accidents to drivers of cars. I discovered that the symbols in the formula actually have 
the following meanings. 
e  A- is the annual accident frequency. 
¢ sand bp are variables whose values depend on the driver's sex: 

$= 0, bp = 3.5 for males; 

= —0.02, by = 2.3 for females. 


 qgand by are variables whose values depend on the driver's socio-economic group: 
g = 0, b; = 13 for drivers in socio-economic groups A, B and C1; 
9 = -0.72, b; = 23 for drivers in socio-economic groups C2, D and BE. 


© pa is the proportion of driving undertaken in the dark. 


© Pb, py and pm are the proportions of driving undertaken in built-up areas, rural 
areas and on motorways, respectively. 


¢  M is the distance driven annually, in miles. 
Ag is the age of the driver, in years, at the mid-point of the accident period. 


X is the number of years since passing the driving test, also determined from the 
mid-point of the accident period. 


I found much of this reassuring, and the rest interesting. The authors of the TRRL 
report don’t talk of chance, but of frequency, of accidents. The idea is that for every 
hundred people (say) of the same age and gender, driving the same number of miles, 
and so on, there will be 100A, accidents in a year. Gender comes into the formula in 
quite a complicated way: it both modifies the overall constant 0.006 33, via the factor 
e*, and has an effect on the way that the accident frequency depends on the number of 
years since passing the driving test. It appears that class (more technically, 
socio-economic group) is a factor, which works in a way similar to gender; I hadn’t 
anticipated this. (Incidentally, the socio-economic groups are defined by occupation, as 
follows: Group A consists of the professional occupations (such as mechanical 
engineer); Group B of intermediate occupations (garage proprietor); Group C1 of 


G. Mayeock, C. R. Lockwood 
and Julia F. Lest 
“The accident liability of car 
drivers’, TRRL Research 
Report 315 (HMSO, 1991). 
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non-manual skilled occupations (driving instructor); Group C2 of manual skilled 
occupations (motor mechanic); Group D of partly skilled occupations (traffic warden); 
and Group E of unskilled occupations (car-park attendant).) 


Another factor I missed is the proportion of driving done in darkness, which on second 
thoughts is clearly important. I find the classification of road types a bit strange, 
though the principle seems clear enough: I think I would have chosen two separate 
classifications, into grade of road (motorway, A class, B class) and situation (urban or 
rural). Note that the variables p,, p- and p» are not independent: since they are 
proportions, and (by implication) one’s driving must always fall into one or other of the 
corresponding categories, they must satisfy py + Py + Pm = 1. 


The TRRL report is rather careful to state the units in which variables are measured, 
where relevant: distances, you will note, are measured in miles. This is an important 
thing to remember when defining variables. 


The time variables Ag and X are also specified very carefully: for example, Ag is not 
just described as the driver's age, but the driver's age ‘in years, at the mid-point of the 
accident period’, It is worth noting in passing that there are limits on the values these 
variables can take: for example, Ag must be a whole number and cannot be less than 
17; there is also an upper limit to its possible values, though this is not so sharply 
determined ~ perhaps 100 would be a reasonably safe choice. It is good practice to look 
out for limitations on values like this. I suppose that X is to be regarded as a measure 
of experience, and to some extent may be considered to stand for skill, which (as T 
remarked earlier) is difficult to measure directly. 


One point about the symbols used: Ag is a single variable, it does not mean A x g. 

I would never use such a symbol myself, because of the risk of confusion. Choosing 
symbols well is a minor but useful art. It helps if each symbol carries some reminder of 
the thing it symbolizes, and it is therefore natural to use the first letter of an 
appropriate word for the symbol. In this case this approach leads to a difficulty because 
two of the main variables, accident frequency and age of driver, begin with the letter 
‘a’, But I think that the use of a or a would have been a better solution to this 
difficulty than the one actually adopted. 


One factor that I raised is omitted from the formula, namely the age and condition of 
the car. I remain puzzled by the absence of any variable relating to the car. This is not 
to say that the TRRL model is wrong; it is as well to be aware from the start that 
mathematical modelling (unlike solving quadratic equations, for example) is not the 
sort of activity in which there is a single correct answer to every problem. In many 
cases there is room for different approaches; this adds to both the interest and the 
difficulty of modelling (in my view, anyway). 


Now that we know what the symbols in the formula mean, we can think about some of 
its implications. Let us look first at the effect of mileage. If all the other variables are 
held fixed, so that we consider the effect of mileage alone, we see that A, is proportional 
to M°279, in other words A, = kM°?" for some k independent of M. Now the graph 
of y = x°279 js similar to that of y = x°*°, or y* = x, which looks like Figure 1. 


Figure 1 
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So, as one would expect, the accident frequency increases with increasing mileage — and 
is 0 if M = 0, which is as it should be! But it is interesting to see just how the accident 
frequency increases. One might have expected that doubling the mileage would double 
the frequency; but the formula predicts that the frequency would be multiplied by 
29279 ~ 1.2, that is, that it would be increased by only about 20%. 


The formula A, = kM°?"® can itself be regarded as a simple model of the frequency of 
road accidents; indeed, it is a perfectly satisfactory model, so long as one is interested 
only in how the frequency depends on mileage. The factor k contains information about 
driving in the dark, the driver's age, and so on; we can tease this information out, at 
the price of making the model progressively more complicated. It is a well-established 
principle of mathematical modelling that it is best to start with as simple a model as 
possible and to make it progressively more complicated by incorporating additional 
features. 


Turning our attention to the effect of driving in the dark, if the proportion of driving 
done in darkness is increased, then the likelihood of an accident increases, all other 
things being equal. If all driving is done in darkness (p, = 1), then the accident 
frequency is 2.6 times greater than if all driving is done in daylight (pa = 0). Ata 
guess, I'd say that in the depths of winter about a third of my driving is done in the 
dark, whereas at the height of summer practically none of it is: according to the 
formula I'm about 50% more likely to have an accident during the winter, which is a 
sobering thought. 


As I had suspected, motorway driving is safer than driving in built-up areas; this 
follows from the fact that the coefficient of p,, is less than the coefficient of Pp; but 
driving on rural roads is safer than both of these. 


Increases in age and in experience both decrease accident frequency, since both Ag and 
X appear in reciprocals, However, the effect of either depends on the gender or class of 
the driver. Why gender is paired with experience, and age with class, is not at all clear 
to me. However, let us look at gender and experience. When all other variables are 
fixed, we have 


Ac = Lexp{s} exp{b2/(X + 2.6)} = Lexp{s + b2/(X +2.6)}, 


where / stands for the product of all the other factors in the formula. For males s = 0 
and bg = 3.5, while for females s = —0.02 and b» = 2.3. It follows that, for a given 
amount of experience, accident frequency is less for a woman than for a man. How 
much less depends on the level of experience. For completely inexperienced drivers, for 
whom X = 0, the ratio of the frequencies (woman to man) is 


exp{—0.02 + 2.3/2.6} : exp{3.5/2.6} ~ 2.37: 3.84 ~ 0.6: 1, 


so the accident frequency for inexperienced female drivers is only about 60% of that for 
inexperienced male drivers of the same age, annual mileage and so on. For drivers with 
20 years* experience, however, the difference is much less marked, the ratio being 


exp{—0.02 + 2.3/22.6} : exp{3.5/22.6} ~ 1.09: 1.17, 
ie. the accident frequency is almost the same for women and men. In fact, the pure 
gender factor e* is relatively insignificant (e~°-°? ~ 0.98), the main difference coming 
from the experience term; as X increases, this difference rapidly becomes less marked, 
as the graphs in Figure 2 show. 
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Figure 2 


I have described several conclusions that can be drawn from the model, and there are of 
course many more. Some of these conclusions merely confirmed expectations, while 
others were in some way new and surprising. The general point to note here is that 
mathematical modelling does not end when one has constructed a model: one has to 
investigate and interpret its implications, both to confirm its validity and to draw new 
conclusions. 


It is time to return to the question of whether all this is a glorious folly. The view 
expressed in The Sunday Times, that it is, presumably rests on two assumptions. One 
of them is stated quite explicitly on page 3 of the magazine: anything that involves 
calculations more complicated than counting up to eight is regarded by the newspaper's 
editors as beyond the capabilities of most of its readers. But there also seems to be an 
implicit assumption that the exercise of trying to analyse the incidence of road 
accidents and the factors that affect their frequency is intrinsically pointless. That 
seems to be the view of the journalists on The Sunday Times; the staff of the Transport 
and Road Research Laboratory presumably thought otherwise. 


In examining the implications of the formula, I have shown that it reproduces some 
well-known or obvious facts (if that is the word) about safety in driving: motorways are 
safer than urban roads, experienced drivers are involved in fewer accidents than are 
inexperienced ones, the more you drive the more likely you are to have an accident. It 
is reassuring that the formula does confirm one’s expectations in this way — if it did 
not, one would have little faith in it (though it is as well to remember that ‘well-known 
facts’ do sometimes turn out to be myths). If the formula was nothing more than a 
complicated way of restating the obvious, then deriving it would indeed have been a 
somewhat pointless exercise. But surely it does contain new information. I was 
surprised by a number of the conclusions drawn above: the effect of darkness, the 
relative safety of rural roads, the way the significance of mileage tails off as mileage 
increases, the fact that experience irons out the difference between the sexes. 


I think that the very personal way in which the formula was presented in The Sunday 
Times (‘a formula to enable you to calculate your mathematical chance of having a 
traffic accident’) both misrepresents it and undervalues it. Though there are some 
pointers to ways in which you and I can reduce our chances of having an accident (try 
to avoid driving in the dark, use motorways and rural roads in preference to urban 


‘These graphs are drawn for 
the same fixed value of 1. 
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ones), they are few in comparison with things over which one has no control. But it is 
scarcely likely that the TRRL intended the formula to be used like this. In fact, the 
formula was obtained as part of a continuing programme of research ‘intended to 
provide a factual basis for more detailed psychological studies’ of the ‘relationship 
between the accident liability of an individual driver ... and factors such as age, sex 
and driving experience’. The TRRL report on this research describes its aims as follows. 


The effect. of age and experience on accident liability are factors of particular 
interest since there are grounds for believing that they influence driving 
behaviour (and hence accidents) differently. Age effects clearly reflect growing 
maturity ~ perhaps changing perceptions of risk, or a growing sense of social 
responsibility. Experience on the other hand reflects a learning process within 
the driving task. The distinction between maturation and learning could have 
important implications for the application of road safety remedial treatments. 
If age is the dominant factor, then short of devising a method of accelerating 
the maturation process, raising the driving age might be the only practical 
countermeasure. If on the other hand, experience plays a crucial part, then 
improved safety might be achieved by devising better ways of imparting those 
skills necessary for safe driving to novice drivers — a matter of training. 


Other possible uses of the results suggest themselves: for example, the analysis of the 
effect of driving in the dark on accident frequency should be useful in the context of the 
debate about summer and winter time, and knowledge of how age affects safety should 
be of interest to insurance companies in determining how premiums should be fixed for 
drivers of different ages. While this may not be the most far-reaching piece of research 
on road safety that has ever been done, to call it a glorious folly does seem to me to be 
a bit unfair. 


It is worth noting that an analysis of this kind itself generates further questions for 
investigation and ideas for improvement. Here, as the authors of the TRRL report 
point out, the nature of the dependence of accident frequency on mileage is crying out 
for explanation. Furthermore, the analysis could be improved by finding a variable to 
stand for ‘experience’ better than the number of years since passing the test. It rarely 
turns out to be the case that there is a final word in mathematical modelling: models 
are almost always capable of extension and improvement. 


1 am conscious that nowhere in this section, so far, have I explained how the TRRL 
formula was actually derived. It was actually done by collecting a lot of data about 
accidents and fitting curves to them, a process related to that described in 

Subsection 1.2 of Unit 3 under the name of empiricism. Here, though, the problem is 
rather more complicated, since there are several variables involved rather than just one, 
and the techniques involved are beyond the scope of this course. The formulation of a 
model is clearly the fundamental step in mathematical modelling, and we must get to 
grips with it eventually; but it is a good idea if first you become more familiar with the 
process of modelling as a whole, and that is the aim of both this and the next section. 


Summary of Section 1 


A mathematical model of a real-world situation or process is an equation, a set of 
equations, or some other mathematical relationship between variables which specify the 
possible states of the system being modelled. In order to understand and analyse a 
model, therefore, it is necessary to know how the variables are defined, and their units. 
Interpreting the model involves working out what it states about the relationships 
between the variables; a successful model will both confirm what one knows already 
about the system under consideration, and predict new conclusions. Models are usually 
developed with a specific purpose in mind, often to solve some particular problem 
relating to the real-world situation being modelled; being aware of this purpose is 
clearly also an important aspect of interpreting the model. 
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2 Two models for you to analyse 


In the previous section I made some general points about the processes involved in 
mathematical modelling, in the context of a description of a specific model concerned 
with the frequency of road accidents. 1 am now going to ask you to confirm your 
understanding of those points by analysing and interpreting accounts of two 
mathematical models. There are two long exercises for you to do in this section: each 
consists of a number of questions about the model described in the preceding report. 

I have selected my examples from prerequisite courses, so you may find one of them 
familiar, but not so familiar that you will not have to think carefully in order to answer 
the questions. 


The first model is taken from M101 and MS284. While you are reading it, bear in 
mind the points made in the previous section, especially the summary. Look out for the 
purpose of the model, the process which is being modelled, the definitions of the 
variables, the derivation of the relation between them, the conclusions that are drawn 
from the model, and any possibilities for improving and developing the model. 


2.1 Modelling drug therapy 


In this report, drug quantities are measured in milligrams and volumes in litres. We Most of this subsection (up 


defined these units when introducing the SI system, but you may like to be reminded to Exercise 1) is taken from 
that M101 Block V Unit 3, 


Section 3.4, or MS284 
1 milligram (mg) = 10~* grams (g) = 107® kilograms (kg) Unit 16, Section 4. 


and that 
1 litre (1) = 10° cubic centimetres (cm*) = 10° cubic metres (m*). 


Patients with asthma suffer from constriction of the airways in the lungs and, 
consequently, difficulty in breathing out. This ailment can be alleviated by introducing 
the drug theophylline into the bloodstream. This is done by injecting another drug, 
aminophylline, which the body quickly converts to theophylline. Once present in the 
blood, however, the drug is steadily excreted from the body via the kidneys. In other 
words, the system ‘leaks’, and unless there is replenishment, the quantity of drug in the 
blood will fall. 


It is known from experiment that theophylline has hardly any therapeutic effect if its 
concentration in the bloodstream is below 5mg/! (milligrams per litre), and that 
concentrations of above 20mg/I are likely to have toxic effects. The problem is to 
administer the drug in such a way that the concentration remains within the 
therapeutic range between 5mg/I and 20mg/I. 


Experimental measurements can be obtained by injecting a known dose of 
aminophylline into a patient, allowing time for it to diffuse throughout the bloodstream 
and subsequently finding the concentration of theophylline in blood samples taken at 
regular intervals. One then knows the initial quantity put into the bloodstream and the 
concentration at various later times, but there is no time for which both the quantity 
and concentration are known. Figure 1 shows a graph of concentration against time 
which could be obtained from such an experiment. 
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Concentration 


(mg/l) 


Figure 1 


The model is based on the following assumptions. 

(a) The conversion from aminophylline to theophylline takes place at the moment of 
injection. This is equivalent to assuming that it is theophylline itself which is 
injected. 

(b) The theophylline is distributed uniformly throughout the bloodstream upon 
injection and at all subsequent times. Hence the drug’s concentration at a given 
time is the same at all points of the bloodstream. 

(c) The quantity of theophylline in the body is proportional to its concentration in the 
blood. The constant of proportionality is called the apparent volume of distribution. 
This is not simply the volume of blood in the body, as the drug is also present in 
other bodily fluids and tissues. It can be thought of as the total volume of blood 
for which the drug concentration would be equal to that in the actual bloodstream. 

(d) No theophylline enters the body between injections. 

(e) The quantity of drug excreted from the body in a given time interval is 
proportional both to the total quantity present at the start of the interval and to 
the duration of the interval. 


On introducing the variables t (time in hours), q (quantity of theophylline present in 
mg) and ¢ (drug concentration in the blood in mg/l), we have q = q(t) and c= c(t). 
Assumption (c) can be expressed as 


a=, 


where V is the apparent volume of distribution in litres. 


Assumptions (d) and (e) are used in concert with the input-output principle, which is 


accumulation =| input | minus | output |: 


The principle is applied to the quantity of theophylline in the body for an interval 
between times ¢ and t + 6t during which no injection is administered. The accumulation 
of the drug over this interval is equal to the difference between the quantities present at 
the end and at the start of the interval, that is, 


accumulation = q(t + dt) — q(t). 
By assumption (d) there is no replenishment of the drug during the interval, and so 
input = 0. 
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Assumption (e) states that the output is proportional both to q(t) and to the duration 
6t of the interval. In fact, the second half of this is only approximately correct, but the 
approximation becomes progressively more accurate for smaller values of ét. Expressing 
assumption (e) in symbols, we have 


output ~ Ag(t)ét, 
where . is a positive constant. The input-output principle, therefore, leads to the 
equation 
a(t + 6t) — q(t) ~ —Aq(t)ét. 
Dividing through by t and taking the limit as ét — 0 gives 
dq 
dt 
This differential equation for g(t) can be supplemented by an initial condition. If a dose 


Dmg is injected at time t = 0 into a patient who was previously free of the drug then, 
according to assumption (a), we obtain 


=. 


q(0) = D. 
Given the differential equation 
dq 
—=- t>0 
et (b= 0) 
and the initial condition 
q(0) = D, 


we must find an expression for q(t), then use the relation 
mi) 
elt) = > 


to deduce an expression for ¢(t). The general solution is 


q= Ae, 


where A is an arbitrary constant, but since q > 0 in the current model we are 
concerned here only with positive values of A. 


Putting ¢ = 0 into the general solution and using the initial condition q(0) = D gives 
D=A, 

so that the required particular solution is 
q=De™. 


Then we have 
D —M 


c=—e™7 


V 
where D/V is the initial concentration. 


A qualitative prediction of the exponential decay solution is that between injections the 
concentration will fall with time but more and more slowly. 


The form of the solution is itself a quantitative prediction, but it can be put into a form 
more amenable to comparison with reality. Taking the logarithms of both sides of the 
equation for ¢ gives 


D 
log, ¢ = log, (?) =e 
This predicts that a graph of log, ¢ against t should be a straight line. 


The qualitative prediction is satisfied by the experimental results shown in Figure 1 
(page 15). The quantitative prediction is also satisfied, for when the experimental data 
are plotted as suggested the result is a straight line (see Figure 2). 
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Figure 2 


To two significant figures, the slope of this line is —0.17 and the intercept on the 
vertical axis is 2.5. Thus we have 


A\=0.17, 


and 


or 


Now D/V is the initial concentration, so for this experiment the initial concentration 
was 12mg/l. For a known experimental dose Dmg, we can calculate that the apparent 
volume of distribution in this case was V litres, where V = 73D. However, this relation 
would normally be used in the form D = 12V to estimate the required dose for a 
patient for whom V is known, when an initial drug concentration of 12mg/1 (in the 
middle of the therapeutic range) is desired, It is found in practice that there is a simple 
empirical relationship between the apparent volume of distribution V (in litres) and the 
patient’s weight W (in kilograms), namely 


V= hw. 
This enables one to estimate the correct dose directly from the patient’s weight. 


For example, for a patient weighing 50kg the apparent volume of distribution is 
3 x 50 = 25 litres, so in order to achieve an ial drug concentration of 12mg/l we 
would need to inject the patient with 12 x 25 = 300 mg of the drug. 


Our model says that after the injection the concentration c is given by 
c= 1267017 


so that it decreases steadily. After about five hours the concentration drops below 
5mg/I (the bottom of the therapeutic range). This could be prevented by administering 
a further injection before this time to send the concentration back up to 12mg/I. 


Such a treatment is feasible but requires many injections (one every four to five hours). 
‘There is an alternative strategy of giving the patient a slow intravenous infusion after 
the initial injection. This continuously replenishes the supply of drug in the blood. The 
mathematical model has to be changed if continuous infusion of theophylline is 
considered. In particular, assumption (d) must be revised, which leads to a new 
expression for the input-output principle and an altered differential equation. 


Exercise 1 


(i) To make sure that you have understood what purpose is served by the modelling 
described in the report, describe in your own words the problem the model is intended to 
solve, how it arises, and why it is significant. Try to concentrate on the essential features 
of the problem, and in describing its significance say whether you think it is specific to 
the situation described in the report, or may have more general application. 


We continue to use the value 
0.17 for \. Experiments show 
little variation in the value of 
A from patient to patient. 
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(iii) 


(iv) 
(v) 
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The problem specified in your answer to the first part of this question cannot be solved 
directly, which is where the model comes in. Fill in the blanks in the following sentence, 
which states what is actually being modelled. 


The model describes how varies with [___] for a given patient, as a 
consequence of excretion via the kidneys. 


Without simply repeating the conclusions of the model, say what form you would 
intuitively expect this variation to take, and why. 


It is assumed in the report that the term ‘concentration’ is self-explanatory. How would 
you explain, to someone who does not regard it as self-explanatory, what it means? 
Your answer should probably contain a reference to one of the ‘assumptions’ numbered 
(a) to (e) which are so conveniently set out near the beginning of the report. 


‘The mathematical formulation of the relationship you have specified in part (ii), for any 
particular patient, involves two variables: define them precisely, and give their symbols. 


What precisely is the model which describes how these variables are related? 


Leaving aside for the moment the question of how the model is derived, let us consider the 
conclusions that are drawn from it. 


(vi) 


Sketch a graph showing how the concentration of the drug varies with time, according to 
the model. Does this accord with your expectations from part (ii)? What evidence is 
provided in the report to show that the model produces realistic results? 


Insofar as the proof of this particular pudding is in the eating, you should be reasonably 
convinced that the differential equation 


dq _ 
Pam: 


gives a satisfactory model of the excretion of the drug. Let us now consider how the model is 
actually derived. The derivation follows directly from a very explicitly stated quantitative 
assumption about the process of elimination of the drug: assumption (e), So everything hangs 
or falls on whether this assumption is a reasonable description of how the kidneys work, 


(vil) Actually, assumption (e) cannot be correct as stated if taken literally, It asserts that the 


amount of drug excreted in a given time interval is proportional to the duration of the 
interval, however long that may be. In deriving the differential equation for q this 
assumption is used only for (infinitesimally) short periods of time. Is the more general 
form of the assumption satisfied by the solution to the differential equation, 

q(t) = De? 


(viii) Assumption (e), even when it is modified, may seem rather arbitrary. To justify it in 


(ix) 


(x) 


(xi) 


detail would involve knowing how the kidneys actually work, and this is a complicated 
business as reference to any physiology textbook will convince you. However, the 
important point is that it is a process of filtration, in which unwanted substances, 
including the drug, are filtered out of the blood. No filtration process is likely to be 
completely effective; but it would be reasonable to assume that every molecule of the 
drug has the same chance of being filtered out during its passage through the kidneys. 
Does this picture of the working of the kidneys lead to the result used in deriving the 
model? 


The model has produced a formula for the variation of the concentration of the drug 
with time, namely ¢(t) = (D/V)e~™, which is both theoretically supported (by the 
argument about the working of the kidneys) and verified by experiment (so far as the 
qualitative predictions are concerned). But can it be used, as it stands, to solve the 
original problem explicitly? What further data, if any, are required to apply the model 
in any particular case, and how are they found? 


The report gives an example of the use of the formula to find when the concentration 
drops below the lower limit of the therapeutic range. According to the report: ‘our 
model says that after the injection the concentration is given by c= 12e~°"" ... after 
about five hours the concentration drops below 5mg/I’. Where does the coefficient 12 
come from, and how is the time of five hours calculated? 


When a patient suffering from a bad asthma attack is admitted to hospital, he or she 
must be put as quickly as possible onto a routine of drug therapy which will keep 
theophylline in the bloodstream at a concentration lying within the therapeutic range, 
say between 12mg/l and 6 mg/l, to be on the safe side. Devise a system of treatment, 
based on the model, which will serve this purpose. What would be the best way to 
present your system of treatment to the nursing staff, given that they have no time for 
complicated calculations in what might be a situation of some urgency, and that the 
relevant characteristics of the patient must be taken into account? 
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(xii) The calculations in the previous two parts have taken an upper limit for the 
concentration of 12 mg/l, while toxicity does not set in until 20mg/I. It is, of course, 
important to err on the safe side, but this looks like excessive caution on the face of it. 
After all, if a larger initial concentration were chosen, the patient could go for longer 
between boosters. Can you think of any reasons why such a low upper limit should be 
taken? 

(xiii) Give two suggestions for further lines of investigation which (if they can be carried 
through successfully) should lead to improvements in this model of intermittent 
administration of a drug. 

(xiv) An alternative to intermittent injection of a drug is continuous intravenous infusion, 
following a single initial injection, This is a method of administering the drug at a slow, 
steady rate. The model discussed above can be adapted to cover intravenous infusion: 
all the assumptions listed still hold, with the exception of assumption (d), which is 
replaced as follows: 

(d) after the initial injection, theophylline enters the body continuously at a constant 
rate. 


Derive (in the form of a differential equation) a model of the administration of a drug by 
intravenous infusion. 


[Solution on page 39 | 


The second exercise is based on an extract from TM282. Unlike the first example, this 
is not designed as a formal report of a piece of modelling, but is more descriptive in 
nature. All the same elements are present, but you may have to work a little harder to 
find them. The basic ideas for this model are covered in MST204 in Unit 4, but in 
slightly different notation: recognizing the appropriate results from that unit is part of 
the exercise! 


2.2 Skid marks 


When the police are trying to reconstruct an accident, the skid marks of a car can be 
very informative. From the length of the skid marks of the wheels, the police can 
estimate the speed at which the car was travelling before the wheels locked and the car 
went into the skid. 


To evaluate the road surface the police sometimes drive a similar ‘test’ car with similar 
tyres, and cause it to skid at the same place, but at a lower speed. They then compare 
the skid marks they produce with the original ones. They assume that the frictional 
forces between surfaces are not dependent upon the speed of the car, only upon the 
mass of the car, the condition of the road and the type of surface. The frictional force 
is assumed to be proportional to the mass of the vehicle, as are any accelerating or 
decelerating forces due to gravity if the car is going up or down a hill. The retardation 
of both cars is the same according to these assumptions. 


So it will be necessary to use the test vehicle's results to calculate the retardation and 
then to use the calculated retardation to estimate the speed of the original car. 


For the test car, I shall call the initial velocity u, and the distance s,. The final velocity 
is 0, so 


uh. 
v=0, 
s= 51 


a is to be found, 


The equation linking u. v, s and a is 


v? =u? + 2as, (+) 
so O=u,? + 2as,, 
au? 
Le. a= 


2s,” 


See TM282 Unit 4, 
Section 5.2. 
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For the original car, if it skidded to rest over a distance s at the above retardation, 


uy? 


s= 


w is to be found. 
Using equation (*) again, 


v? =u? + 2as, 


s 
So u=h 
ETS 


So the speed of the car before the accident can be estimated. Note that the driver's 
original speed would be greater than this if the car had not stopped at the end of the 
skid. 


Exercise 2 
(i) State in your own words the purpose of this model, and say when it may be useful. 


(ii) What is the role of the test car? Why could the police not just be issued with tables 
giving the speed in terms of the length of the skid? 


(iii) The description of modelling that I have been trying to put across in this unit is that its 
basic concern is with finding a relationship between variables that specify the system 
under consideration; in other words, to find a formula which will enable you to calculate 
something if you know the value of something else, or which tells you how something 
varies with something else. What are the appropriate ‘something’ and ‘something else’ in 
this model (in words)? The second sentence of the description of the model tells you! 


On the basis of what you know about skids, say what you can, in the simplest terms, 
about the nature of this variation. 


(iv) The term ‘retardation’ is used several times in the description of the model, without 
explanation. What do you think it means? Is it a positive or a negative quantity? 


(v) There are quite a few symbols floating around in the description of the model, but 
according to part (iii) we are really concerned with the relationship between just two 
variables. It seems to me that the symbols could conveniently be grouped into three; 
symbols for the fundamental variables (of which there are two); symbols for the data; 
and symbols which are introduced to make the calculations easier. Classify the symbols 
according to this scheme. 


(vi) No units of measurement are given anywhere in the definitions of the symbols. This is 
not the practice adopted in MST204, and you are advised not to follow it in this course. 
But does it matter in this instance? 


(vii) What is the basic model which underlies the whole discussion? 


(viii) Unlike the report of the drug therapy model, the assumptions that lead to the model are 
not explicitly listed here, though they are described in the text. It seems to me that 
there are really two fundamental assumptions: one is concerned with a special property 
of the acceleration (or retardation) of a car which is skidding, which permits the use of 
the formula in the answer to the previous part of the exercise; the other says something 
about the relationship between the crashed car and the test car, which enables one to 
use the measurements of the test car to make deductions about the crash. Can you 
formulate two such assumptions? 


(ix) The assumptions actually made in the text support, or give details of, these two basic 
assumptions. State what support there is in the text for the assumption of constant 
acceleration, and what support there is for the assumption that the crashed car and the 
test car have the same acceleration. 


(x) How do Newton’s second law and the properties of kinetic friction given in Unit 15 
justify assumptions made in the text (for a level road)? (You might find it helpful to 
draw a force diagram for the skidding car — indeed, I hope that this idea has occurred to 
you without my hint.) 
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(xi) The possibility that the crash took place on a slope is alluded to in the description of the 
model, but it is not stated explicitly whether the model applies to such a situation. Does 
it? You should be able to answer this question by considering whether that situation is 
covered by the assumptions. On the same basis, decide whether the model would apply 
to the following situations: a crashed Rolls Royce (assuming the police cannot afford 
such large and expensive cars); a crash into a very strong headwind; a crash in a shower 
of rain, if the road dries out before the police arrive on the scene. 


(xii) The model, so far as the crashed car is concerned, might be summarized as follows. 
We want to determine the initial speed u in terms of the length of the skid s (which is 
the distance the car travels before coming to rest). Under the assumption of constant 
acceleration this is given by 

ue+2as=0 or u=V/—2as. 
(The car’s acceleration, a, is negative.) 
Confirm that this result agrees with your intuitive expectations as given in part (iii), 
Why is this not the end of the story, and where does the test car come in? 


(xiii) In the application of the constant acceleration formula to the original car, the equation 
= appears. This does not look very informative: what do you suppose it means? 


(xiv) By combining the formulae for the crashed car and the test car, the final formula 


8 
u=u 
it 


is obtained. Says the text: ‘so the speed of the car before the accident can be estimated’. 
If the actual skid marks are four times as long as those of the test car, and the test car 
was going at 40 mph when it skidded, was the crashed car breaking a 70mph speed limit? 


(xv) Suppose that you are a police instructor: explain to a policeman about to go out on an 
accident investigation for the first time how to estimate the speed of the car before the 
accident. 


(xvi) Explain the final sentence: ‘the driver's original speed would be greater than this if his 
car had not stopped at the end of the skid’. 


[Solution on page 42] 


Summary of Section 2 


Reading and interpreting descriptions of models can be hard work: there is usually a 
lot of information to absorb, and it can be difficult to focus on what is really important. 
For example, it seems to me that in the drug therapy model, the name of the drug does 
not help one’s understanding of the mathematics (though it is of crucial importance to 
the clinician, of course). I often think of the process of coming to terms with a model as 
one of excavating: digging away until you find what lies at the root of it all. In the drug 
therapy model it is the differential equation —Aq which lies at the root, whereas in 
the skid marks model it is the rule relating initial and final velocity with distance, for 
motion with constant acceleration. The great skill in modelling (as you will see when 
You come to construct models for yourself) is to formulate the fundamental equation or 
relationship which describes the process in which you are interested. There are several 
pointers in these two reports to how this is done. First, there is the importance of 
picking out appropriate variables and defining them carefully. Secondly, it is necessary 
to simplify matters in order to make progress: for example, recognizing in the drug 
therapy model that the mixing problem could be postponed, if not ignored. Thirdly, it 
is good practice to record the assumptions one makes in deriving the model: it certainly 
helps the reader, but more importantly, it clarifies things for you and it suggests ways 
of developing the model later. Fourthly, note the importance of the collection of 
relevant data, both to check the predictions of the model and to furnish the values of 
any constants that are needed to apply the model to any particular situation. 
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3 The skills of modelling 


One thing should be clear to you by now, if it wasn’t already clear before you started: 
mathematical modelling involves many different skills. To be good at mathematical 
modelling you have to be able to do some mathematics, of course; but you have to be 
able to do many other things as well. If you tend to think of mathematics as a set of 
procedures (such as the procedures for solving differential equations) then you may not 
think of these additional skills as being relevant to mathematics. However, most 
practising mathematicians, concerned with pure mathematics as well as applied, would 
regard that interpretation of mathematics as a travesty of the truth. One of the most 
powerful motivations for studying mathematics is the desire to solve new problems, in 
other words, problems for which there is no known solution procedure. The skills 
required in mathematical modelling are instances of general problem-solving skills. 

To be able to solve mathematical modelling problems is more useful, and more difficult, 
than being able to solve first-order differential equations by the integrating factor 
method. It calls on skills of invention, analysis and interpretation which apply to all 
sorts of problems, not just mathematical ones (pure or applied). 


Here is a list of skills that may be required in the solution of a modelling problem. 
Interpreting the question. 

Choosing appropriate variables. 

Simplifying the problem and identifying its most important features. 
Specifying relationships. 

Interpreting and testing the results. 


Improving the model. 


When one is tackling a modelling problem in earnest, one has to call on these skills 
repeatedly, in a complicated and interactive way. In particular, ‘improving the model’ 
will usually involve repeating the whole modelling process, and therefore using all the 
other skills over again. It seems wise to practise them individually first; this is the 

aim of this section. Each subsection deals with one of the skills from the list. The 
individual skills are illustrated with examples taken from the models in the last section. 
These illustrations are followed by exercises for you to do, most of which are based on 
models developed in previous units of this course. 


3.1 Interpreting the question 


In mathematical modelling the problems one faces are rarely posed in a way which can 
be translated directly into mathematical form. For example, from the description of the 
treatment of asthma in the report of the drug therapy model, one might have thought 
that in order to construct a model one would have to know quite a lot about drugs. 
However, we were able to produce a straightforward statement of the underlying 
problem which could serve as the basis for a model, namely 


to find how the quantity of a drug in the body varies with time. 


It is important to bear in mind thé purpose of the model. The purpose of the drug 
therapy model is to devise a timetable for the administration of the drug which will 
keep its concentration between the therapeutic limits with the least disturbance of the 
patient. 


This is typical of the approach that is required in order to get started on a modelling 
problem. 
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Exercise 1 

The effects of double glazing were investigated in Unit 12. Figure 1 shows the temperature 
variation through a double-glazed window. 

Based on this diagram, and the methods of heat transfer which produce such a temperature 


variation, a formula for the rate of loss of heat through the window can be derived: it takes 
the form 


9 = AU (Bin — Bout), 


where q is the rate of heat loss in watts, A is the area of the window in m?, U is a constant 
called the U value of the window, with units Wm7?°C~', and @jq and Gont are the 
temperatures of the air inside and outside the building, in °C. 


Summarize the problem which is investigated in the double glazing model, and suggest why 
one might want to carry out such an investigation. 


Exercise 2 


In Unit 15 we considered the athletics field event of putting the shot. By analysing the motion 
of a shot, considered as a particle moving in a vertical plane under the action of gravity alone, 
the range R of the shot is found to satisfy the equation 
g 2 m 

Ta? couto® + (tan@)R+h=0, 
where vp is the speed with which the shot is launched, @ is the angle at which it is launched, 
measured from the horizontal, and h is the height of the point from which it is launched above 
the ground. 


What contribution could mathematical modelling make to the improvement of a shot-putter's 
performance, and how? 


Exercise 3 

Many shops do good trade in selling sandwiches for office workers’ lunches. But no one likes to 
buy stale sandwiches, and quite a few fillings become unsafe to eat in a short time. The 
manager of a sandwich shop clearly has a problem when it comes to knowing how many 
sandwiches to provide for sale each day. Can you specify the problem in a form which is 
sufficiently precise to serve as the starting point of a mathematical model? Can you think of 


any other commodities, besides foodstuffs, whose sale leads to a similar problem? 
[Solutions on page 44 | 


3.2 Choosing appropriate variabl 


As has already been said, the identification of the main variables is of paramount 
importance. If you can summarize the problem in terms of determining how one 
quantity varies with another, or several others, you should have no difficulty in 
identifying the variables. In the case of the drug therapy model there are two main 
variables: the first is a measure of the quantity of drug (either the mass of drug in the 
body, or the concentration of drug in the bloodstream); the other is the time since the 
drug was administered. In the skids model the problem is to find how the length of the 
skid depends on the speed of the car when skidding began, from which the variables are 
obvious. It pays to be careful in defining variables (‘time since the drug was 
administered’, not just ‘time’). One should give units, even though I haven’t done so 
here, 


Exercise 4 
In Unit 4 we posed the following problem. 


‘An object which is initially at rest is dropped from the Clifton Suspension Bridge and falls 
into the River Avon, 77.0 metres below. Assuming that the only force acting on the object is 
the force of gravity, find the time of fall.” 


Select two main variables for this problem, and define them carefully. 


Unit 12, Subsection 4.3 


Bin 


Figure 1 


Unit 15, Section 4 


Unit 4, Subsection 3.3 
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Exercise 5 


‘The modelling of a mortgage was discussed in the first TV programme of the course. The 
programme dealt with the repayment of a loan of £20000 over a period of 25 years at an 
interest rate of 15% per annum. A recurrence relation for u,, the amount owed after r years, 
was obtained: 


Urs = 1.15u, — 12M, 
where M is the monthly repayment, 
What are the main variables in the mortgage problem? 


Exercise 6 


New supermarkets open so frequently that I can't help wondering how the managements of 
supermarket chains decide where to site new stores. One fundamental factor in the choice of 
site must be the catchment area. Can you define variables which would be useful in setting up 
a mathematical model which could be used to decide the position of new supermarkets? 


[Solutions on page 44 | 


3.3. Simplifying the problem and identifying its most important 
features 


The skid model depended on the assumptions that the acceleration of a skidding car is 
constant, and that for given conditions different cars have the same acceleration while 
skidding. These assumptions follow from a well-established theory of sliding friction, 
but hold only if (for example) air resistance is ignored. To ignore air resistance is 
justified on two counts: first, its effects are probably small compared with those of 
sliding friction; secondly, the resulting model is relatively easy to analyse. In modelling 
one should always look for as simple a model as possible consistent with the salient 
features of the problem (to have ignored the effects of friction completely would 
obviously have been counter-productive). But it is important to be clear about the 
simplifying assumptions that have been made in order to arrive at that model. This is 
necessary to ensure that anyone who may be interested can follow how the model has 
been derived, Furthermore, if you have an explicit list of the assumptions you have 
made, and you wish to improve your model, you have an obvious place to start: review 
the assumptions and ask which can and should be modified or relaxed. 


Exercise 7 


What simplifying assumptions would you make in setting up the equation of motion in the 
Clifton Suspension Bridge problem described in Exercise 4? 


Exercise 8 
What assumptions must be made to derive the equation q = AU(@, — one) in the double 
glazing model of Exercise 1? 


Exercise 9 


How much of a car's windscreen or rear window is actually cleaned by the windscreen wiper or 
wipers depends on the length of the wiper blade(s) and the positioning of the wiper(s) on the 
windscreen. It is obviously desirable that the area cleaned should be as large as possible. 
What simplifications would you make in order to set up a model of the operation of 
windscreen wiper(s) that would help to determine the best arrangement of the wiper(s)? 


[Solutions on page 45 | 


3.4 Specifying relationships 


The use of the input-output principle to derive the fundamental differential equation in 
the drug therapy model is a good illustration of quite a common technique; we used the 
same technique in Unit 3 when developing models for populations. Of course, most 
models of mechanical systems involve the use of Newton’s laws, though the use may be 
hidden, as in the skid model. 


Unit. 1, Subsection 4.1 
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Exercise 10 


Explain, in your own words, how the recurrence relation in the mortgage model (described in 
Exercise 5) is derived. In the problem considered in Unit 1, specific values for the interest rate 
and the loan are given, But all repayment mortgages work in the same way, so it should be 
possible to write down a recurrence relation for any such mortgage. Can you do so? You will 
have to introduce into the model some parameters which you will have to define. 


Exercise 11 


In Unit 8 a highly simplified model of the suspension of a car was discussed. The basic 
features are shown in Figure 2. The car body is represented by a particle, which is supported 
by a perfect spring and a damper, together representing the suspension. To test the 
suspension, a potential buyer of the car depresses the body a certain distance below its normal 
position (with the car stationary) and releases it. 


Describe how the equation which determines the subsequent motion of the car body is 
obtained. 


Exercise 12 


A supermarket where I frequently go shopping is surrounded by car parks as shown in the plan 
below. 


NZ 


When I go there I can never decide where to park. The problem arises because the closer to 
the centre one goes, the longer one has to spend driving around looking for an empty parking 
space; if one parks far away from the centre, however, one has to spend a long time walking to 
the shops. It seems to me that I could try to resolve this difficulty as follows. I shall try to find 
how the time between arriving at the car park and walking into the shopping area depends on 
where I park, and then work out how far away to park to minimize this time. Let ¢ be the 
time between arriving at the car park and walking into the shopping area, and x the distance I 
park from the shops. The time it takes to find somewhere to park is a decreasing function of 2: 
I suspect it is inversely proportional to x. 


Formulate a relationship between t and x that could be used to solve the problem. 


Figure 3 


[Solutions on page 45] 


Unit 8, Subsection 1.4, 
Exercise 7 


car body 


damper 


Figure 2 
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3.5 Interpreting and testing the results 


Here there are several points to be emphasized. It is usually worth checking that the 
model predicts the kind of results that one would expect from common sense and from 
experience. If possible, one should validate the model by comparing its predictions with 
data from an experiment or other reliable source. It is good practice to try to 
reformulate the results so that this check turns into something simple like drawing a 
straight line, as in the drug therapy model. The model may not directly provide the 
information you want to solve the initial problem, so you may have to do some further 
mathematics (solve a differential equation, or find a maximum or minimum, for 
example). But you may also require some data to give an explicit numerical solution to 
the problem ~ such as the value of \ in the drug therapy problem. The beauty of the 
skid problem lies in the use of the test car to provide these extra data. 


Exercise 13 


Unit 3 dealt with the growth of populations. An initial attempt to model the growth of a 
population, based on the principle that the change in the population in a short interval of time 
is proportional to the population size at that time, leads to the equation 

dP 

oa AP, 
where P is the population size, t is the time, and \ (a constant) is the proportionate rate of 
increase of the population. To obtain a more sophisticated model, it is assumed that the 
proportionate rate of increase of the population is not constant, but decreases linearly as the 
population size increases. ‘The resulting equation for the population size is the logistic 


equation, 
dP P 
> (1-3) 
(where a and M are constants), which has the solution 


a M 

~ 1+(M/Po = 1en*t* 

with Pp the initial population. Data for the growth of a population of yeast cells, when 
plotted, produced the following graph. 


P 


Population 
700 


Time (hours) 


Figure 4 


What qualitative features of the solution to the logistic equation suggest that it is an 
appropriate model for population growth? How can the solution be rearranged to provide a 
straightforward way of checking the model against data such as those given for the population 
of yeast cells? 


Exercise 14 


The recurrence relation for the mortgage problem quoted in Exercise 5 was used to find the 
value of M for which the loan would be repaid exactly after 25 years; the value obtained was 


M = 257.8324... 


Unit 3, Section 2 
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Tt was discovered that a monthly repayment of £258 (obtained by rounding the value of M to 
the nearest pound) leads to an eventual overpayment of £427. 


Describe exactly how the recurrence relation for the mortgage problem is used to find the 
monthly repayment, taking into account the practical problems of implementing the solution 
(taking the effect of rounding as an example) as well as the mathematical ones of deriving it. 


Exercise 15 


An attempt to model cooking times, which was the subject of a television programme in M101 
and MS284, leads to the following result. Consider a spherical body of mass m, initially at 
room temperature, which is placed in an oven maintained at a much higher temperature. 
The aim of the model is to find the time it takes for the centre of the body to reach a 
predetermined temperature (the temperature at which food is deemed to be cooked). 
The model predicts that this time t is given by 

t=km™, 
where k is a constant. What checks can you devise to confirm that this is a reasonable result? 
How would you use the result to find haw long it would take to bake a potato? 


[Solutions on page 46] 


3.6 Improving the model 


The points that were ignored in simplifying the model can now be reconsidered. In the 
case of the drug therapy model it was suggested that the initial mixing of the drug with 
the bloodstream should be considered. A less demanding modification of the skid 
model was briefly raised in the final part of the exercise, which asked about the effect of 
relaxing the implicit assumption that the car stops at the end of the skid. Since the 
purpose of the exercise is to find incontrovertible evidence that the driver was speeding 
when the accident occurred, a qualitative answer was enough. 


Exercise 16 


Explain why the exponential model of population growth, described in Exercise 13, is usually 
inadequate. Discuss the modifications which lead to the logistic model. 


Exercise 17 


In the Clifton Suspension Bridge problem (Exercise 4), an initial model is proposed in which 
velocity-dependent forces are ignored; this model is superseded by a better one in which such 
forces are taken into account. What is the effect of this on the result? 


Exercise 18 


When a drug is taken by mouth, rather than being given by injection, there is a delay before it 
gets into the bloodstream because it has to pass to the stomach and be absorbed, Can you 
devise a model of the change in concentration over time of a drug taken by mouth? 


[Solutions on page 47] 


Summary of Section 3 


The solution of a modelling problem will usually require the use of the following skills. 
e Interpreting the question. 

¢ Choosing appropriate variables. 

e Simplifying the problem and identifying its most important features. 
e Specifying relationships. 

e Interpreting and testing the results. 

e Improving the model. 


Identifying and practising these skills at every opportunity will pay dividends when you 
come to do some modelling for yourself. 


M101 Block V Unit 1, 
Section 1.4 
MS284 Unit 14, Section 4 
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4 Putting it all together 


In previous sections you have learnt how to interpret models and you have practised the 
individual skills required in modelling. You must now begin to integrate these skills, so 
that you can tackle modelling problems yourself. In this section I am going to talk my 
way through a modelling problem, describing how I would think about it; I hope that 
this will give you some guidance as to what is involved in doing one’s own modelling. 
Most of this section consists of reading only, though there are some exercises for you to 
do at the end to help you consolidate the points that have been made; but I hope you 
will read in an enquiring frame of mind, and keep asking yourself at each stage what 
you would do next if you were tackling the problem. My general line of attack is similar 
to that already brought out in the discussions of the drug therapy and skidding models 
in Section 2; so you should probably find yourself always one step ahead of me as you 
read this passage. 


In the next section you are asked to put the ideas into practice, by doing a piece of 
modelling yourself, with guidance. 


The question I am going to consider is this: 
Immediately after a cup of tea is poured out, it is usually too hot to drink. 
How long does it have to be left to cool before it is drinkable? 


This is not, on the face of it, a terribly important question, and falls within the 
category of questions one asks out of curiosity rather than because the answer is likely 
to be useful. However, there are many different situations in which cooling plays a 
significant role, some of which are of considerable importance; you should not find it 
difficult to think of industrial processes in which cooling plays an important part, for 
instance. The cooling of a cup of tea is a familiar and simple example, which involves 
the principles which apply in most cooling problems. 


lam now going to describe the steps I would go through in trying to answer this 
question, 


‘The problem is concerned with how hot a cup of tea is: the first thing to do is to find a 
precise way of stating what this means. How does one measure how hot something is? 
With a thermometer. In other words, the problem is concerned with the temperature of 
the tea. Furthermore, the question asks how long it takes for something to happen, so 
time is going to be of major importance. The guts of the problem, therefore, is to find 
how the temperature of the tea varies with time, and the two main variables can be 
identified already: they are 


4, the temperature of the tea, in °C; 
t, the time since the tea was poured out, in seconds, 


My interpretation of the problem is that @ is a function of f, and I have to find out as 
much as I can about this function — if possible, I want to express @(t) explicitly in terms 
of t. Limagine that, if I can do so, the expression will contain various constants (which 
I tend to call ‘parameters’) which represent various relevant physical characteristics of 
the system: the initial temperature of the tea and the temperature of the surroundings 
are obvious examples. More will no doubt present themselves as my ideas develop. 


Before I try to analyse the cooling process in detail, I shall ask myself how, intuitively, I 
would expect the temperature to vary with time. The first and most obvious point is 
that anything which is hotter than its surroundings cools down: so the temperature of 
the tea will progressively decrease as time passes. How can I express this 
mathematically? By saying that 6(t) will be a decreasing function. But can the 
temperature keep on decreasing indefinitely? It is surely a matter of experience that 
however long one leaves it, one’s tea never gets colder than its surroundings; but if one 
waits long enough, its temperature will get as near that of the surroundings as makes 
no difference. What would a graph of temperature against time be likely to look like? 
Figure 1 shows a couple of possibilities. 
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Figure 1 


Which is the more likely? The second is more likely than the first, surely: physical 
processes tend to be smooth. 


It occurs to me that I could make an educated guess at the function represented by the 
second graph. Its shape is that of a decaying exponential, but it is moved up a little, 
parallel to the vertical axis. Can I use this to guess a formula for @(t)? It suggests that 
the relation between @ and t takes the form 


A(t) =a+be™*, 


where a, b and are constants. Here be~™* is the basic exponential: to obtain the 
shape shown in the figure, b and A must both be positive. The addition of the constant 
a corresponds to moving the exponential up parallel to the vertical axis. This formula 
for @(t) gives a partial answer to the question ‘how does @ depend on t?’; it would be of 
more use, however, if I could interpret a, b and A in terms of the parameters of the 
problem. This can easily be done for a and b by thinking about initial and (if you see 
what I mean) final conditions. First, putting t = 0, I obtain 


(0) =a+b. 


But (0) is the initial value of the temperature of the tea, its temperature when it is 
poured out; it would fit in with the usual practice if I were to denote this by 0). On the 
other hand, I suggested that in the long term (as t + oo) the temperature of the tea 
will tend to that of its surroundings, say @,; but as t 00, e~™’ — 0, so this tells me 
that a = @,. I can now see that b = 0) — @,, and so 


A(t) = 0, + (00 — O,)e-™*. 


So far so good, but \ is still undetermined. What is more, the formula was obtained by 
a piece of heroic guesswork, How could I confirm the accuracy of my guess, and how 
can I find the value of A? One way would be by experiment: take a cup of tea and 
measure its temperature periodically. If I were to plot the results (measured 
temperature against time) I would hope to find that the points lay on a graph like that 
in the previous figure. How could I check the model more effectively? By making a plot 
of log,,(0(t) — @,) against t, because the model predicts that the result would be a 
straight line, of slope —A. Not only would this give confirmation of the model, but it 
would provide the value of 4. I expect you will have noticed the similarities between 
this and the drug therapy model. 


This procedure has its limitations, however. It would confirm the result for a particular 
cup of tea, and might reasonably be supposed to provide supporting evidence of the 
correctness of the model for cups of tea in general. But because it is basically empirical, 
it is of rather limited validity. My aim must surely be to find a method of predicting 
the time of cooling for any of the many situations in which the knowledge might be 
useful, and having to carry out a new set of measurements in each case cannot really be 
described as a satisfactory approach. 


It would be better to find an analytical derivation if possible. To do so I must think 
more carefully about the process of cooling. So at this point I must make a mental 
effort and try to recall what there is in Unit 12 that might be relevant and useful. That 
unit is about the flow of heat, and I haven't thought about heat at all yet in my 
consideration of the problem. The first question is ‘how are heat and temperature 
related?’. The essential points seem to be these: where there is a difference of 
temperature, heat will flow (from the warmer to the cooler place); and when something 
loses heat, and the loss is not made good, then its temperature will fall. So the 
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following picture of the processes underlying the cooling of the tea begins to form in my 
mind: since the temperature of the tea is higher than that of its surroundings, it loses 
heat continuously; since the tea is losing heat, its temperature falls continuously. The 
greater the difference in temperature between the tea and its surroundings, the greater 
the rate at which heat is lost, and the greater the rate at which the temperature of the 
tea falls. Conversely, as the tea cools, it loses heat more slowly, and the rate at which it 
cools slows down. In the long term, there is practically no difference between the 
temperature of the tea and that of its surroundings, no loss of heat, and no further 
cooling. It sounds as though this mechanism will reproduce the qualitative features of 
the cooling process that have been discussed already, which is promising. So next I 
shall try to turn these ideas into a mathematical model. 


Three ways in which heat flows from one place to another are discussed in Unit 12, as I 
am sure you will remember: they are radiation, conduction and convection. Which of 
these will apply in our case? Radiation is connected, in my mind at least, with things 
like electric fires and the Sun; if I put my hand near a cup of tea, I won't feel any 
warmth worth speaking of unless I actually touch the cup or the tea, so there can’t be 
much heat being radiated and it seems safe to ignore it. Both of the other methods of 
heat transfer will apply, though: heat is conducted through the cup, and is convected 
directly from the surface of the tea and also from the outer surface of the cup. To keep 
things simple it would be a good idea to concentrate on just one of these forms of heat 
loss (conduction or convection), if possible: this will be justifiable if it seems likely that 
heat is lost much more rapidly by one route than the other. The effect of ignoring the 
less significant of the two will be that the answer may be a little inaccurate, but it is 
better to develop a simple model which you can be confident will actually produce an 
answer, than to develop one that is so complicated that you cannot use it. 
Improvements can be made bit by bit, but only if you have something which actually 
works to improve on. 


The outside of a cup containing hot tea is a lot cooler than the tea itself: after all, one 
reason why china, for example, is used to make cups is that it is a good insulator. So 
heat is probably lost much more rapidly by convection from the surface of the tea than 
it is by conduction through the cup. To begin with, therefore, I shall ignore all forms of 
heat loss except convection from the surface of the tea, Now, what do I need to know 
about the tea and its surroundings, in addition to the difference in their temperatures, 
in order to calculate the rate of loss of heat by convection? I need to know the area of 
the surface of the tea: let me call this Am?. 1 also need to know the convective heat 
transfer coefficient: let me call this hWm~* °C". Then if q is the rate of loss of heat 
(in watts), and the temperature difference between the tea and the surrounding air is 
4-9,, 


q=hA(O-6,). 


That has made a start on the first half of the cooling process: I have a model of the 
way the temperature difference between the tea and its surroundings leads to a flow of 
heat from the tea to the surroundings. The other half of the process is concerned with 
the resulting fall in the temperature of the tea. Near the beginning of Unit 12 there is 
a formula which relates a change in the internal energy of a body to the consequent 
change in its temperature. You may remember this formula, which is given in the form 


Ey = Ey = me(02— 61). 


In Unit 12 this formula is discussed in the context of the heating of water in a kettle: 
the heat transferred to the water by the kettle raises its internal energy from E to E2, 
with resulting increase in temperature from 0; to 42; m is the mass of the water, and c 
a constant, characteristic of the substance being heated, called its specific heat. This 
formula applies equally to the fall in temperature resulting from a loss of energy. In our 
case, it can be rewritten as follows: when a mass m kg of tea loses 5E joules of heat 
energy, its temperature falls by 6E//(mc) °C, where ¢ is the specific heat of tea, in 
Tkg- C71. 


I now have to put the two parts of the process together. There is, however, a problem: 
the formula for the rate of loss of heat by convection applies only in the steady state, 
that is, when the temperatures concerned remain constant in time; this is definitely not 


Strictly speaking, this 
formula holds only in the 
steady state. We shall deal 
with this point shortly. 


Unit 12, Subsection 1.3 
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true of the situation under consideration — indeed, the whole point of the exercise is to 
find how the temperature of the tea varies with time. But I can get round this by a 
standard device, which you have met once already in this unit. I shall consider what 
happens in a short interval of time of length ét, during which the temperature of the 
tea remains effectively constant. I can then use the steady-state result as an 
approximation. It is an approximation which improves as 6t — 0. The price I have to 
pay is that I shall finish up with a differential equation. 


In the interval of time from t to t + ét, the rate of loss of heat by convection, using the 
steady-state approximation, is 
a(t) = hA(O(t) — 4,). 
The quantity of heat lost is therefore 
5E = q(t)dt = hA(O(t) — 6,)ét. 
When the tea loses this amount of heat energy its temperature drops by 6E/(mc), and 


so. 


A(t + 6t) = a(t) — = = 0(t) — ao) ~0,)6t. 


It follows that 
A(t + 6t)— A(t) A 


3 a (04t) ~0,). 
Now let 6t — 0 to obtain 

do hA 

Fal). 


This is my analytically derived model of the cooling process. 


Once you have a differential equation, there is one obvious thing to do: try to solve it. 
But before doing so, one might ask whether it predicts the right kind of behaviour for 
@. We expect. 4 to decrease with t so long as @ > @,, and to decrease more slowly the 
nearer @ is to @,. Is the differential equation consistent with this behaviour? Well, so 
long as 0 > 0,, we have d@/dt < 0, so @ is decreasing; and the closer that @ is to 0,, the 
smaller is |d0/dt|. Note that @ = @, is a constant solution, or equilibrium state (like the 
terminal velocity for a particle falling under air resistance, or the equilibrium 
population size for a population growing according to the logistic equation). So the 
equation does appear to predict the right kind of behaviour. 


How is it solved? One way is to use the integrating factor method (the integrating 
factor is e™*, where \ = hA/(me)). Another is to use separation of variables (divide 
through by @—0,). Here is a third method, quicker than either of the others. Put 
@=0-0,: then dé/dt = d0/dt (because 0, is constant), so the equation becomes 

Be en = Ng, eines 

dt me me 
an equation whose solution you should know (if only because we have met it before in 
this unit, in the drug therapy model): 


Ot) = de™. 
In terms of the original variable @ we have 
hA 
O(t) —0, = (0) —O,)e™™"*, = A= — 


I should just like to be sure this makes sense. The formula agrees with the one I 
guessed before by thinking about the graph of @ against t. The temperature @ tends to 
9, in the long term. The temperature will fall more quickly if the surface area A is 
increased: that seems right, because in that case there will be more room for the heat 
to escape. On the other hand, if you have more tea it will cool more slowly: a greater 
mass will hold more heat, and the loss of a given amount of heat will result in a smaller 
drop in temperature. So the answer does appear to make sense. 


It remains to answer the original question: how long will it be before the tea is 
drinkable? How do we set about answering this? So far as purely mathematical 


These ‘equations’ are, 
strictly speaking, only 
approximations; however, 
they get more and more 
nearly true as ét gets smaller. 


If I had anticipated this 
outcome, I would have used 
@ as my main variable in the 
first place! 
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processes are concerned, it is simply a matter of setting 6(t) equal to the temperature 
at which tea becomes drinkable, say @,, and solving for t: this gives 


04 — 4 = (Oo — A.)e-™, 


Oa —Os _ ye 
Ga. e 
s tot ogy (Oe) sie (ae 
e = 88 \G= as) RA NBO 


But this isn’t the end of the story, is it: in order to obtain an explicit numerical answer, 
the values of the various parameters in this formula need to be known. Here a list 
would be useful — or, even better, a table like the following. As well as listing the 
meanings of the symbols and the units of measurement, I have asked myself in each 
case how the value is to be obtained (as no doubt you have too). 


Symbol Definition Units How obtained 
m Mass of tea kg Weigh a full and an empty cup 
—use cold tea to avoid accidents! 
c Specific heat of tea Jkg~!°C>! Use value for water given in Unit 12 
h Convective heat transfer coefficient Wm7?°C~' Look up in book of physical 
between tea (water) and air constants 
A Area of surface of tea m? Measure diameter of top of cup 
A Initial temperature of tea BO Measure with thermometer: 
presumably just below boiling point 
64 Temperature of surrounding air bi 2 Measure with thermometer 
O4 Temperature at which tea becomes Lc Sample, then measure with 
drinkable thermometer 


Once the values of these parameters have been found, the rest is arithmetic — or is it? 
Here are some real-world values (actually for a mug rather than a cup). 


Quantity Value 


m 0.25 kg 
c 4200 Jkg=! C=! 
h 250 W m=? "C=! 


A 0,005 m? 
% 80°C 
0, 20°C 
04 60°C 


With these values my calculator gives the following answer: 
t = 340.590 6908. 


It is very tempting to conclude that the time taken for the tea to cool enough to be 
drinkable is 340.590 690 8 seconds, or 5.676511513 minutes. What would be wrong in 
doing so? I have derived a model, based on good physical laws and principles; I have 
correctly calculated the result this-model predicts, using reliable data; so the answer 
must clearly be right. But can one really trust an answer given to 10 significant figures 
— an answer that purports to have a far greater degree of accuracy than any watch on 
the market? The data are given to at most two significant figures, so we can scarcely 
expect the result to be more accurate. In any case, the model itself is approximate, as 
we know: the effects of conduction, for example, have been ignored. It would be much 
more honest to state the result as ‘the time taken for the tea to cool enough to be 
drinkable is about 6 minutes’. 


The consequence of ignoring the effects of conduction should also be considered. The 
values of convective heat transfer coefficients are difficult to determine accurately, and 
they depend on the precise conditions in a way that is difficult to predict: this is 
implicit in the enormous ranges of values quoted in Table 3 on page 17 of Unit 12. 
Whether or not the tea is stirred, and whether or not there is a draught in the room in 
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which you are drinking it, could seriously affect the value of h. I suspect that the 
correction resulting from taking account of conduction through the cup would be less. 
than the error resulting from the indeterminacy in the value of h, so that from the 
practical point of view there would be little to be gained from developing the model in 
that direction. However, it is worth asking whether the inclusion of the effects of 
conduction would change the qualitative behaviour of the model, and whether it would 
tend to increase or decrease the estimate of the cooling time. Since the formula for the 
rate of heat loss by conduction takes the same form as that for the rate of heat loss by 
convection ~ in each case the rate of heat loss is proportional to the difference between 
the temperatures of the tea and the surrounding air — the only change to the model will 
be that A will be replaced by another term. Moreover, the new term will be larger 
than /A (since it will incorporate additional heat loss), so the effect of the inclusion of 
conduction in the model will be to decrease the estimate of the cooling time slightly, as 
one would expect, 


I do not seriously suggest that whenever one has a cup of tea one should calculate how 
long to leave it before it is cool enough to drink. The benefits of the foregoing exercise 
were intended to be entirely pedagogical. You should have learnt two things from it. 
One is a new (so far as this course is concerned) piece of information about heat, 
namely the law that governs the cooling of hot bodies. This is known as Newton’s law 
of cooling, and can be stated briefly as follows: the rate of decrease of temperature is 
proportional to the temperature difference. The other is some further insight into the 
nature of modelling. I'd like to summarize the points that have arisen which, in my 
opinion, throw useful light on the nature of modelling, in the form of a collection of 
questions that one is likely to have to ask oneself at various stages when one is engaged 
on a modelling problem. 


Questions to ask yourself while modelling 


What process will I have to model in order to answer the question posed? 


Here the question asked was ‘how long does it take for a cup of tea to cool?’; the 
process to be modelled is the transfer of heat between an object and its surroundings, 
and the consequent changes in the object’s temperature. I put it in this rather vague 
way so that you can see that the model is actually much more general than the 
particular situation apparently under investigation: the same principles should apply to 
the warming of cold objects as to the cooling of hot ones, for example. 


What are the key variables I shall have to use in my model? 


“How does one measure how hot something is?’ turned out to be a useful pointer to 
temperature, one of the two key variables in this case. It is very easy to overlook time 
as a variable, simply because it is so obvious: it is almost a tautology to say that time is 
likely to be a key variable in a problem concerned with how long it takes for something 
to happen. Likewise, spatial variables (coordinates, for example) are bound to be 
required in problems concerned with spatial variation or change in position. In most of 
the modelling problems you have met or will meet in this course, the key variables will 
be few in number — usually two or three. (The accident frequency model with which 
this unit started is untypical in this respect — but it was chosen to make a point.) 


Can I formulate the underlying modelling problem in terms of how one of 
my variables depends on the other(s)? 


If you have a clear idea of the key variables you will probably have already answered 
this question, and vice versa; the only extra point is to decide which is the dependent 
and which the independent variable, as we say in the trade. My succinct formulation of 
the cooling problem is ‘how does the temperature of the tea vary with time?’; here 
temperature is the dependent variable, time the independent one. Forms of relationship 
other than functional dependence arise in modelling problems, of course, but the 
question, suitably modified, is still relevant. 
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What do common sense and experience tell me about the qualitative 
aspects of the relationship between my key variables? 


I got quite a long way with the cooling problem by this form of guesswork. You won't 
always be so lucky, but thinking about the problem in this manner may suggest, for 
example, a hypothesis to test empirically, which will be some sort of progress even if 
you can’t see how to develop a theoretical model. In many cases (such as the accident 
frequency problem) this will be the only way of getting anywhere. When you can 
develop a theoretical model, on the other hand, observations along the lines of *y must 
be an increasing function of x’ provide useful qualitative checks. (But always remember 
that neither common sense nor experience is completely reliable; indeed, most of the 
really important advances in science have been due to someone showing that what 
everyone else thought was common sense actually isn’t.) 


What would I learn from an experiment? 


If you are contemplating some kind of empirical approach, you must be clear what you 
are looking for. Experiment, or the collection of masses of data, for its own sake, is not 
good modelling. There are usually two reasons for experimenting, or collecting data by 
some other means: one is to test a hypothesis about the nature of a relationship 
between variables; the other is to find the values of parameters which occur in the 
model. In the cooling problem an experiment could have been used to confirm that the 
temperature falls exponentially with time, and also to find the value of the parameter \ 
which measures the rate of cooling. These are examples of the two different possible 
uses. 


How can I find a relationship between my variables theoretically? 


There is no simple answer to this question, no way of producing a Procedure Box for 

modelling. But there are several tips and suggestions that I can give. 

e First, it can often help if you have some idea from the start what kind of model you 
are expecting to obtain. An ordinary algebraic equation? An equation involving 
trigonometrical functions, or exponentials? A differential equation? If so, what will 
its order be? 


e Look for relevant laws or principles, such as the input-output principle, the law 
governing heat transfer by convection, Newton's laws, and so on, The idea of 
considering ‘what happens in a short period of time of length ét and letting 6t tend 
to zero’ should be mentioned here too. 

e Ask yourself whether you have previously met other modelling problems with 
similar features (the drug therapy problem has surprisingly many similarities with 
the cooling problem, for example). 

Try to develop a simple picture of the processes involved, in your own terms and 
language. This could take the form of a metaphor (the kidneys are like a filter). A 
diagram is often useful — essential, I would say, in the case of mechanics problems, 
to show the forces. 

« Simplify like mad. In most real-life situations there are many processes going on 
simultaneously. Some will be quite irrelevant to the problem in question (whether 
the beverage is tea or coffee); others can be ignored, at least initially, because their 
effects are relatively insignificant (radiation and conduction compared with 
convection). Your aim should be to identify just the one or two most important 
mechanisms at work, and model those. If this leads to a somewhat superficial and 
inaccurate model, this doesn’t matter very much, since it is always possible to try 
and incorporate additional features later. It is much better to have a simple model 
which you can subsequently try to improve, than to be prevented from making a 
model at all because you are trying to cover everything. 


© Lists can be helpful, such as the list of the different mechanisms for heat transfer 
mentioned in the solution to the cooling problem. You can work through the items 
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on the list and decide which have the most significant effects. If you can put them 
in order of importance, so much the better (convection first; conduction second; 
radiation nowhere). I have to admit that this order would be more convincing if I 
had done some rough calculations to estimate the relative sizes of the rates of loss 
of heat by each of the three methods. But, assuming that I'm right, a clear policy 
emerges: model the heat loss by convection first, incorporate conduction if 
necessary next, but don’t bother with radiation at all. There is another benefit to 
be gained from drawing up such a list: when reporting on your modelling to 
someone else it is very helpful to the audience to be explicit about which 
mechanisms you have taken account of, which you have ignored, and why — in other 
words, to make explicit your assumptions in the manner of the report on the drug 
therapy model ~ and a list will help you to do this. 


« Don't be afraid to introduce new variables, such as q in the cooling problem, and 
parameters such as h and A. (The distinction here is between the variables @ and t, 
which would appear in any cooling problem, and the parameters, h, A, m, c, ,, 
which are constants for any particular cooling problem, but whose values vary 
between different cooling problems. The rate of loss of heat, q, and the energy 
change, 5E, are variables, but have a different status from @ and t: q and 6E are 
introduced to make the derivation easier, but they are hidden, or subsidiary, in the 
sense that they don't appear in the final model.) As you think about the processes 
involved, the need for these new variables and parameters should become apparent. 
Don’t let them multiply unchecked, however: if you are running out of letters of the 
alphabet, you probably haven't simplified enough. Do keep a list of symbols you 
introduce, with their definitions and units: there is nothing more irritating than 
deriving a beautiful, complicated formula and then realizing you have forgotten 
what a, means. 


When I have derived a model, what do I do next? 


Now is the moment to remind yourself what the problem is that you are trying to solve, 
and to reformulate the problem in terms of the mathematical model you have obtained. 


How do I use my model to solve the original problem? 


Do some mathematics! If you have obtained a differential equation, or a recurrence 
relation, solve it. If your problem asks you to find a maximum or minimum, 
differentiate and set the derivative to zero. Of course, what mathematics you do will be 
determined by the problem that was originally posed. Remember that numerical 
answers are usually required in the real world, and so this part of the exercise involves 
thinking about how one finds the values of any parameters that are needed to obtain a 
numerical answer. Remember also to give the answer to a level of accuracy which is 
appropriate to the problem, the model and the data, to the best of your judgement; in 
particular, don’t give the appearance of claiming greater accuracy than the model and 
the data can support. 


How can I be confident I’m on the right track? 


This is a good point at which to make some qualitative checks. Confirm that your 
model predicts the kind of qualitative behaviour you expect, including, for example, 
what happens in the long term and in limiting cases (as I did in the air resistance 
problem and the logistic model). Try the effect of varying the value of a parameter, and 
seeing if the answer makes sense (as I did for m and A on page 31 in the cooling 
problem). 


How can I check that the model corresponds to reality? 


You might consider conducting an experiment, as described in the drug therapy model 
for example; or you might look for confirmatory data from elsewhere, such as reports of 
other investigations of a similar problem. 
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How can I improve my model? 


Here is where your list of assumptions will be useful: look at what you decided to leave 
out the first time round, and consider how it can be incorporated in the model. 


In Unit 3 we used the following diagram to summarize the various activities involved in M101 Block V Unit 1, 


modelling. Section 1.4 
MS284 Unit 14, 
Subsection 4.1 
Figure 2 


‘The list of q: estions and answers above constitutes an attempt to describe in greater 
detail what . ay be involved in the first six of these seven activities. (We shall discuss 
writing reports in the Guide to Mathematical Modelling.) The diagram can help in 
modelling, but is intended only to remind you of the things that you may have to do; as 
we pointed out in Unit 3, it is not intended to be a straightjacket. 


Here are a couple of problems on the theme of heat transfer, for you to think your way 
through to be sure you have understood the points which have been raised in this 
section, You are not being asked to solve the problem in either case, just to think how 
you would set about solving it. (For the sake of completeness, the solutions to these 
exercises contain more detail than you are expected to provide in answering them.) 


Exercise 1 

1am fond of fresh milk, provided that it’s cool: I do not gain much pleasure from drinking 
milk near room temperature, I keep my milk in the refrigerator, of course: I want to know for 
how long a bottle of milk remains pleasurably drinkable after it has been taken from the 
refrigerator. How would you set about solving this problem? 


Exercise 2 

A thermostatically controlled domestic central heating system works roughly as follows, When 
the temperature in the house falls to a predetermined lower value, the heating is turned on, 
and delivers heat at a constant rate, until the temperature rises above a predetermined higher 
value. The heating then goes off and the honse cools, until the cycle begins again. How would 
you set about modelling this situation? 


[Solutions on page 49 | 
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Here is a list of points it may be helpful to consider when trying to develop a 
mathematical model. It is based on the questions that I have suggested you ask 
yourself while modelling, but I have grouped them so that they correspond to the 
stages of the modelling diagram. 


1. Specify the real problem 
* Decide what process you will have to model in order to answer the question 
posed. 
2, Set up a model 
e Pick out the key variables that will have to be used in the model. 
e Try to formulate the underlying modelling problem as thow does one of my 
variables depend on the other(s)?’. 
e Use your common sense and experience to say what you can about the 
qualitative aspects of the relationship between your key variables. 


e Formulate a relationship between your variables, by using relevant physical 
laws, principles like the input-output principle, your own mental picture of the 
processes involved, and any similar models you have met before, remembering 
always to simplify the problem as much as you can. 


¢ State what simplifying assumptions you have made in formulating your model. 


3. Formulate the mathematical problem 
e  Re-interpret your initial problem in terms of your model; convert the modelling 
problem into a mathematical one ... 
4. Solve the mathematical problem 
e ... and solve it. 


5. Interpret the solution 


¢ Decide what data you will need in order to specify the values of any 
parameters, and collect them. 


e Make any qualitative checks you can think of, to confirm that you are on the 
right track. 


¢ Interpret your mathematical solution so as to give a solution to the original 
problem in appropriate terms. 
6. Compare with reality 
e Decide what data you will need to validate your model, and collect them. 


e = Think of ways of improving your model, by relaxing one or more of your 
simplifying assumptions. 
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5 A modelling exercise for you to tackle 


What shape should a tin can be? ‘Cylindrical, of course,’ I expect you say (but is it 
really ‘of course’?). But given that tin cans are cylindrical, what shape of cylinder: 
short and fat, long and thin, or somewhere in between? I’m thinking of the standard 
kind of can in which soup and beans and other foodstuffs are preserved. Most of the 
ones in my kitchen cupboard contain around 400 grams of food and are rather taller 
than they are wide. There's a degree of uniformity in the general shape of cans which I 
find quite surprising, when there are so many different brands: even imported cans 
seem to be the same shape as home-manufactured ones. Is this tradition, or is it 
because they all conform to some ideal shape; and if the shape is not ideal, could there 
be advantages in changing it? 

In this section you are asked to investigate the best shape for a tin can. This is a 
modelling problem for you to do, but you are not left entirely on your own to tackle it: 
the problem has been broken down into steps. 


Exercise 1 
Consider the best shape for a tin can to contain a given amount of baked beans. 


(i) Tt all depends what you mean by ‘best’ .... What do you think the word ‘best’ should 
mean, in the phrase ‘the best shape for a tin can’? 


(ii) ‘Try to formulate a clear statement of a modelling problem out of this. (A clear 
statement doesn't contain imprecise words like ‘amount’ !) 


(iii) Specify your variables. 


(iv) How does the volume constraint affect things? 


If you don’t worry too much about seams and joins, you should find it fairly straightforward to 
write down a formula which will help you to solve the problem, 


(v) Do so, expressing the area A in terms of radius r and volume V. 


(vi) List the assumptions that you made in deriving this formula. 


It is worth pausing and checking that this formula makes sense. 


(vii) What happens to A as r — 0, and as r — 00? Interpret your answers in terms of the 
shape of the can. 


(viii) Find the shape of a can of minimum area, Summarize your answer in words. Are tin 
cans really this shape? 


(ix) What have we left out? 


(x) Modify the formula for the area of tin-plate used to make a can so as to take into 
account your answer to part (ix). Solve the modified problem. Comment on your answer. 


(xi) Would square section cans use less tin-plate than is used in the optimal solution found in 
the previous part? 


[Solution on page 51 | 
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Appendix: Solutions to the exercises 


Solutions to the exercises in Section 2 


1. (i) The problem is to find how frequently it is 
necessary to give successive doses of a drug used in the 
treatment of a disease. The problem arises because the 
concentration of the drug in the patient's bloodstream 
decreases over time as the drug is excreted via the kidneys; 
the drug ceases to be effective if its concentration falls below 
a lower limit, but there is an upper limit to the size of a dose 
because the drug is toxic if its concentration is too high. The 
problem is important because the administration of the drug 
involves injection, which is painful, and so it is desirable to 
find how to administer it so as to ensure safety and 
effectiveness while avoiding unnecessarily frequent injections. 
To have a rule for this, in terms of relevant characteristics of 
the patient, would be better than simply waiting until the 
patient became ill again before repeating the injection. 

It seems to me that, in this answer, I have been able to 
describe the essential features of the problem without 
naming the disease or the drug, and without having to 
discuss the details of the way it is administered. I think, 
therefore, that the specific problem considered in the report 
is just one of a general class concerned with the use of drugs 
in the treatment of a wide variety of diseases, and indeed 
with other situations such as the administration of 
painkillers. May there not even be situations entirely 
unconnected with medicine where similar processes are at 
work: what about the chlorination of water in a swimming 
pool? 

(ii) The model describes how the quantity of drug in the 
body varies with time, for a given patient, as a consequence 
of excretion via the kidneys. 

Or: 


The model describes how the concentration of drug in the 
bloodstream varies with time, for a given patient, as a 
consequence of excretion via the kidneys. 

The model does actually predict, in addition, how the 
quantity (or concentration) of drug at a fixed time depends 
on the patient's apparent volume of distribution, and on the 
size of the initial dose. However, we are primarily concerned 
with the time, since the question we have to answer is how 
frequently to repeat the dose; in the interests of simplicity 
we should concentrate on the time dependence first (and 
when we do so, the other factors enter naturally anyway). 
The quantity of drug (and its concentration) will clearly 
decrease with time; moreover, the less drug there is, the less 
concentrated it will be, the less effective will the kidneys be 
in removing it, and so the more slowly it will decrease. 


(iii) It is assumed that the drug ‘is distributed uniformly 
throughout the bloodstream ... at all ... times’ 
(assumption (b)). This means that, at any given time, the 
mass of drug contained in any sample of blood is directly 
proportional to the volume of the sample; the concentration 
is the constant of proportionality. To put it another way, the 
concentration of the drug is the mass of drug in a sample of 
blood divided by the volume of the sample; the assumption 
of uniformity says that this will not depend on the size of 
the sample nor on whereabouts in the bloodstream it is 
taken, but only on the time at which it is taken. ‘Sample’ 


here must be taken in a hypothetical sense: I am describing 
what would happen if one did take a sample, rather than 
supposing that an actual sample is taken. 

Concentration is analogous to density. 

In a more complicated (and more realistic) model one might 
relax this assumption, and consider what happens if the 
drug is distributed unevenly through the bloodstream. It 
would become necessary to define concentration in a more 
sophisticated way. Concentration will depend on position in 
the circulatory system. It is defined as the mass of drug in a 
sample of blood taken from the neighbourhood of the given 
position, in the limit as the volume of the sample tends to 
zero. 

When the drug is first administered it is concentrated close 
to the point at which it enters the bloodstream; as the blood 
circulates it carries the drug along wit! mixing it with 
the blood until it becomes effectively uniformly distributed. 
In the model developed here, that process is ignored, or if 
you prefer, regarded as taking place instantaneously. Note 
how fair the report is in pointing this out as an assumption, 
rather than leaving the reader to work it out. 


(iv) There are two possible answers to this question, 
depending on how you answered part. (ii). If you chose the 
quantity of drug, then your answer to this question should be 
q the quantity — better, mass — of drug in the body, 
in milligrams, 
t the time since the drug was administered, in hours. 
If you chose the concentration, on the other hand, then your 
answer should be 
ce the concentration of drug in the bloodstream, 
in milligrams per litre, 
t the time since the drug was administered, in hours. 


The relationship between q and c requires a little thought. 
The equation e({t) = q(t)/V appears to be a restatement of 
the definition of concentration given above (under the 
assumption of uniform distribution), but is not quite, 
because g(t) is defined as the quantity of drug in the body, 
not the bloodstream, and as pointed out in assumption (c), 
the drug is distributed to parts of the body other than just 
the bloodstream. Both variables, q and c, are required: q 
because its initial value is the size of the dose, and we can't 
tell how this is divided between the bloodstream and the 
rest of the body; ¢ because it is the concentration which 
determines the effectiveness (and the toxicity) of the drug. 


Note also that the specification of the time origin as the 
moment of administration of the drug depends on the 
assumptions that the conversion of aminophylline to 
theophylline, and the distribution of theophylline uniformly 
through the bloodstream, take place instantaneously 
(assumptions (a) and (b)). 

As long as we are discussing a single patient, the quantity V 
(the apparent volume of distribution) should be considered 
as a constant rather than a variable; and as long as we are 
discussing what happens to a single dose, the size of the 
initial dose D should also be considered as a constant. 
Quantities like V and D, which are not fundamental 
variables, but whose values determine the conditions in any 
single instance of the model, are often called parameters. 


40 


(v) For my money. the best answer to this question is the 
differential equation relating q and t: 


This is, after all, a direct mathematical description of the 
process of excretion of the drug via the kidneys. I suppose a 
case could be made for describing the formula for the 
concentration, 

Dow 


as a model, but the term is usually reserved for the most 
basic relationship between the variables arising directly from 
the assumptions, rather than those which may subsequently 
be derived by applying mathematical techniques, such as 
solving a differential equation. 

The differential equation could have been reformulated in 
terms of c before any attempt was made to solve it, instead 
of solving it first and then substituting for q in terms of c, 
This would have given a model in terms of ¢ and t: 


where D, V and . are positive constants. The graph of ¢ 
against ¢ is thus a negative exponential, and looks like this. 


Figure 1 


‘The evidence in the report that this is the correct formula 
for the concentration is based on the results of a (somewhat 
idealized sounding) experiment to measure the concentration 
of drug in blood samples taken at regular intervals. (Note 
that the description of the experiment states that time is 
allowed for the drug to diffuse throughout the bloodstream: 
we should not expect the experiment and the predictions of 
the model to match initially, while the real drug is not 
uniformly distributed, because this is not accounted for in 
the model.) If ¢ depends exponentially on ¢, then the plot of 
Jog, ¢ against t should be a straight line, and the 
experimental results confirm this, 


(vii) The question reduces to consideration of whether the 
equation 

a(t + 6t) — q(t) = —Aq(t)dt, 
derived from the assumption by use of the input-output 


principle, holds exactly when q(t) = De~™, when 6t is not 
necessarily small. Well, in this case 


g(t + 6t) — g(t) = De"** — De = De“ (e™™* — 1), 


whereas 
—q(t)it = —ADe~™*6t. 

Thus 
q(t + 6t) — g(t) A —Aa(t)dt. 

The assumption is modified in the course of finding the 


output, where it is pointed out that the proportionality of 
the output to the duration of the interval is only 


MST204 M Solutions 


approximately correct, though the approximation improves 
as 6t decreases. (In the final report of this model, to be 
found in the Guide to Mathematical Modelling, the 
assumption will be corrected to reflect this modification.) 


(viii) Yes, it does. Let the rate of flow of blood through the 
kidneys (which I assume to be constant) be p litres per hour, 
Then in a time interval of length t hours, the volume of 
blood flowing through the kidneys is pét litres, and the mass 
of drug passing into the kidneys is thus approximately 
e(t)pét mg. (Note that this deduction holds only for short 
periods of time, so that the variation of e(t) with time t can 
be ignored: hence the use of the word ‘approximately’.) 

The number of drug molecules in this quantity of drug is 
proportional to the mass of drug, and since each molecule 
has the same chance of being filtered out, a fixed proportion 
of this mass of drug is excreted. Thus 


q(t + 6t) — q(t) = —ke(t)pét, 
where k is the constant of proportionality, But 
e(t) = q(t)/V, from which it follows that 


a(t + 5t) ~ q(t) = ~Agq(t)ét, 
with A= —kp/V. 


(ix) In order to apply the formula for the concentration to 
any particular patient, one needs to know the values of D, V 
and A. There is no difficulty about D, which is the dose: it 
is under the control of the doctor. But V and A cannot be 
measured directly: to find V directly one would first have to 
drain off the patient's blood and measure its volume; and 
even that process, impractical as it is, would not provide the 
answer since V is the apparent volume of distribution, 
Hence the importance of the indirect method, in which the 
values of V and . are deduced from the plot of log, ¢(t) 
against ¢. This is a straight line, with slope — and intercept 
log,(D/V) on the vertical axis, as follows from the equation 


log, (2) — At. 


Since the value of D is known, the value of V can be 
calculated. 


Even so, this is not a really satisfactory solution so far as 
clinical use is concerned. A seriously ill patient can't be 
subjected to the kind of experiment described in the report 
before the drug which will ameliorate his or her condition 
can be administered. There are two facts, additional to the 
model as discussed so far, which make it useful in drug 
therapy; both are mentioned explicitly in the report. First, 
there happens to be a simple relationship between apparent 
volume of distribution and body weight, so that the value of 
V can be deduced in a straightforward way from the result 
of weighing the patient. Secondly, the value of \ does not 
vary very much from patient to patient. Thus if we write 


log. ¢ 


c(t) = Been, v=iw, 


we have a pair of formulae which can be used to find when 
the concentration drops below the therapeutic level and 
which requires only a knowledge of the size of dose and the 
patient's weight for its application. 


(x) The coefficient 12 is simply the initial concentration of 
the drug (since it is e(0), the value of c when t = 0). It is 
tempting, but unnecessary, to derive this by calculating D 
and V, In fact, the relationship between ¢(0), D, V and W 
is much more useful as a method for calculating what dose is 
required to give a desired initial concentration (which must 
be kept below the toxicity level), rather than the other way 
round. 
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As t increases from 0, c(t) decreases from 12. To find when 
the value of ¢ drops to 5 we must solve the equation 

Bx 12e7Ot 
for t. The solution is 

p= 1oBe(12/5) _ 0.88 

O27” (0:17 

so the time for the concentration to fall to the lower limit of 
the therapeutic range is indeed about 5 hours. 


(xi) There are three steps. 


=5, 


1. Find what dose is required to produce an initial 
concentration of the drug of 12mg/I. 

2. Find how long it will be before the concentration falls to 
6 mg/l, when a further injection will be necessary. 

3. Find what dose is required to bring the concentration 
up to 12mg/I again. (A smaller dose will be required at 
this stage because even though the concentration has 
fallen to the bottom of the therapeutic range, there will 
still be a considerable amount of drug in the body.) 


Further injections of the dose found in step 3 may then be 
repeated at the intervals found in step 2 for as long as the 
treatment is required. 
‘The initial dose Dmg is given by 

D/V=12, where V= iW, 
W being the patient's weight (in kg), so D = 6W. 
To find the interval between further injections, we proceed 
in the same way as in the previous part, except that now we 
have to solve the equation 


UES beds 
‘The solution is 
= Mee? 2.00", 
*= oar o1r 


Thus after about 4 hours the concentration will have 
dropped to about 6mg/I. To boost it up to 12mg/] again 
requires just half of the initial dose, that is, 3W mg. 

The instructions to the nursing staff might look as follows. 
1. Weigh the patient on admission. 


2. Give an initial dose of 60mg for every 10kg of the 
patient's weight. 


3. Give a booster dose of half the initial dose every four 
hours. 


This assumes that there is a sufficient safety margin in 
working well within the therapeutic range. as we have done, 
that we need consider the patient's weight only to the 
nearest 10 kilograms. It is certainly a good idea not to 
complicate nursing procedures unduly, by demanding that 
the dose be calculated correct to ten places of decimals, or 
that the period between injections be timed to the nearest 
second. Nursing procedures must be pretty well-conditioned 
to small changes in the data (to adapt the terminology of 
numerical analysis) to be practical, as anyone who has been 
in hospital will know. 


(xii) The answers to previous parts of the question point 
out some of them. There is bound to be a fair degree of 
uncertainty in any model of a bodily function, and in 
particular the value of the constant will vary somewhat 
from patient to patient: someone who eliminates the drug 
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relatively slowly is in danger of overdosing if the upper limit 
is set too high. Moreover, the relation between apparent 
volume of distribution and body weight is not a precise law. 
It is unrealistic to expect that the doses will be repeated at 
exactly four-hourly intervals, however well run the hospital 
is. Indeed, the therapeutic limits themselves may differ 
between individuals. 

However, there is one consideration which probably 
outweighs all these, and it relates to the assumption that the 
drug is instantaneously mixed uniformly through the 
bloodstream. All the concentrations calculated in the model 
depend on this assumption. Now there is clearly a period of 
time, immediately after injection, when the drug is not 
uniformly distributed but is concentrated in a small region 
of the bloodstream; and so initially the true concentration, 
in this region, is much higher than the calculated one. It 
follows that quite a considerable degree of erring on the side 
of caution is necessary if toxicity is to be avoided. 


(xiii) The two which are suggested by the answers to the 
previous two parts of this question are as follows. 


1. Carry out an ‘error analysis’, to find out how sensitive 
the conclusions drawn from the model are to variations 
in the value of A, the relation between V and W, the 
measurement of W, and the time between injections. 


2. Model the mixing process which takes place after the 
injection while the drug is being distributed through the 
bloodstream. 


You may have thought of others. 
I do not pretend that either of my suggestions is easy to 
carry out! 
(xiv) The model will consist of a differential equation for 
the quantity of drug in the body as a function of time, as for 
the intermittent case. A new symbol must be introduced, to 
stand for the rate of infusion: suppose it to be J mg/hour. 
The input-output principle may be applied; the version 
given in the report still applies to the excretion of the drug 
by the kidneys, but now the drug is being continuously 
replenished. In a time interval of length ét hours, an amount 
of drug [ét mg is administered, i.e. 

input = /ét. 
Thus the input-output principle gives 

a(t + dt) — q(t) ~ Tt — Aq(t)ét, 
whence 
dq 
B= F= dg. 
F alad 
Note that this differential equation has a constant particular 
solution q = I/A; ideally, one would like to choose I so that 
the value of J/) corresponds to a concentration which lies in 
the middle of the therapeutic range, since a suitable initial 
dose will then ensure that the concentration remains always 
at this level. In fact, the general solution of the differential 
equation is 

g(t) =1/A+ke™ 
for some constant k, from which it follows that q(t) tends to 
1/2 in the long run; so this is a very good regime because it 
tends to settle down to the desired behaviour of its own 
accord, and the initial dose doesn’t have to be very accurate, 
as long as it is below the toxic level! 
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2. (i) | When a car skids, its tyres may leave marks on the 
road. The purpose of the model is to find a method of using 
the lengths of such skid marks to work out how fast the car 
was going when it began to skid. This information may be 
useful to the police after a road accident, when they are 
trying to find out what happened. (It would be interesting 
to know whether an estimate of the car’s speed obtained like 
this would be acceptable in a court of law as reliable 
evidence in a prosecution for dangerous driving or a similar 
offence, To answer this question, both legal and 
mathematical knowledge would be required. A 
mathematical modeller would not be expected to have the 
necessary legal expertise; but he or she should be prepared 
to tackle the mathematical issues, of which the most 
important are the reliability and accuracy of the results.) 


(ii) The test car provides data which are needed for the 
calculation, in addition to the length of the skid marks in 
the accident. These data allow a direct comparison to be 
made between the length of the skid marks in the accident, 
for which the speed of the car is unknown, and the length of 
skid marks made by a car whose speed at the onset of the 
skid is known. In effect, the data for the test car are used to 
work out the deceleration of the car involved in the accident 
while it is skidding. This is possible because the conditions 
under which the test takes place are similar to those for the 
accident, apart from the speed of the car. The state of the 
road surface and the car’s tyres, and other circumstances of 
the accident, will have a significant effect on the skid. They 
will vary a lot from accident to accident, It is easier to 
reproduce the conditions in a test than it is to make 
allowance for them in a table. 


(iii) The relationship required is one which gives the speed 
of the vehicle at the onset of the skid in terms of the length 
of the skid marks. The longer the skid marks, the faster the 
car must have been going when it began to skid (all other 
things being equal), so the speed will be an increasing 
function of the length of the skid marks. 


(iv) It sounds as though ‘retardation’ should mean the 
same as ‘deceleration’. However, in the middle of the 
passage the phrase ‘the above retardation’ is used to refer to 
the formula a = —u,?/(2s,), so that retardation is equated 
with acceleration. It is therefore negative in this case, 
(Much confusion can be avoided by sticking to the term 
‘acceleration’, and remembering that for a point moving in 
the positive direction its acceleration is positive when its 
speed is increasing, negative when its speed is decreasing.) 


(v) The fundamental variables are the original car's initial 
speed, u, and the length of the skid, s. The data, which are 
the corresponding quantities for the test car, are represented 
by the symbols u (initial speed) and s; (length of skid). All 
the symbols appearing in the final formula have now been 
identified, but there are still others, which occur only in the 
course of the calculations. These are a, the acceleration, and 
v the final speed, of either car. The final speed is zero in 
each case (so the symbol v is not really necessary). The 
reasons why the two cars may be assumed to have the same 
acceleration will be considered in parts (viii) and (ix). The 
question of units (or rather, the lack of them) is dealt with 
in the next part of the exercise. 


(vi) No, it doesn’t matter in this instance. One way of 
seeing why involves rewriting the final formula as 
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‘The left-hand side of this equation is the ratio of two speeds, 
while the right-hand side is the square root of the ratio of 
two distances. Now the ratio of two speeds is simply a 
number, and has the same value regardless of whether the 
speeds are measured in ms~', mph or leagues per century 
(provided that both speeds are measured in the same units). 
Likewise, the ratio of two distances, takes the same value 
whatever unit of measurement is used. The omission of units 
in the description of the variables is therefore not important 
in this case. Some units must be used when one is making a 
measurement, of course; but provided the same units are 
used for the speeds, and the same units for the distances, it 
does not matter which. One could even use different systems 
of units for the two kinds of measurements: the formula will 
give the speed of the original car in mph if the speed of the 
test car is expressed in mph, even if the distances are 
measured in metres, However, mixing units is a practice 
which should be avoided unless you are really sure you know 
what you are doing. 

(As you will know if you have taken TM282, the omission of 
units is not an oversight. In that course the approach to 
units of measurement is different from that taken in 
MST204, and in some ways is more sophisticated. It is 
based on the fact that any useful and significant formula 
takes the same form in any system of units, provided that 
consistent units are used throughout, Thus the formula 
relating initial and final speed, acceleration and distance for 
a point moving with constant acceleration is the same 
whether SI or imperial units are used, We therefore do not 
need to specify the units for the variables v, u, a and s in 
the formula uv? = u* + 2as, provided we use consistent units 
for distance, speed and acceleration.) 


(vii) The model which underlies the whole discussion is the 
standard model of the motion of a point moving in a 
straight line with constant acceleration, which you met in 
Unit 4, Section 1. The formula required here is the one 
relating the initial and final speed, acceleration and 
distance. In the box on page 15 of Unit 4 this is given as 


v? = va + 2a0(x — x0), 


which becomes the same as the equation used here when w is 
substituted for vp, a for ag and s for x — x. 


(viii) The basic assumptions are as follows. 
1. A skidding car has constant acceleration, 


2. The test car has the same acceleration as the original 
car. 


(ix) The support I would give for the two assumptions is as 
follows. 


1. The text says that *...the frictional forces between 
surfaces are not dependent upon the speed of the car, 
only upon the mass of the car, the condition of the road 
and the type of surface’. This implies that the frictional 
force is constant. If the road is horizontal, it follows 
immediately that the acceleration will also be constant 
(assuming that friction is the only force acting). If the 
road slopes at a constant angle, the same will be true, 
since the additional component of force due to gravity 
will also be constant. 


2. The next sentence says that the ‘frictional force is 
assumed to be proportional to the mass of the vehicle, 
as are any ... forces due to gravity if the car is going up 
or down a hill’. The fact that the test car has similar 
tyres to the original car, and the conditions of the test 
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being similar to those of the crash, mean that the other 
factors affecting the frictional force are the same in the 
two instances. The total force component acting is 
therefore a constant times the mass of the car, where 
the constant is the same for the test car and the original 
car. Since acceleration is force divided by mass, this 
justifies the assumption that the two cars have the same 
acceleration. 


(x) The rather wordy solution to the previous part of this 
exercise does suggest that a diagram will clarify what is 
going on. For an accident which takes place on a horizontal 
road, assuming that air resistance can be ignored, the force 
diagram looks like Figure 2. 


7 
if 
E 


Here N denotes the magnitude of the normal reaction and F 
the magnitude of the frictional force. By Newton's second 
law, 

ma=—F, 


—> direction of motion 


Figure 2 


where m is the car's mass and a its acceleration in the 
direction of motion, Furthermore, since there is no vertical 
motion, 

N=mg. 
Since the car is skidding (with its wheels locked) the rule 
relating frictional and normal reactions, given in the box on 
page 30 of Unit 15, applies. This gives 


F=,'N, 
where ji’ is the coefficient of kinetic friction. So 
a=—p'g. 


This says that the acceleration of a skidding car is constant, 
and that the magnitude of the acceleration is independent of 
the mass of the car and will be the same under all conditions 
leading to the same value of the coefficient of kinetic 
friction. The test is designed with the precise aim of 
producing conditions leading to the same value of the 
coefficient of kinetic friction. 


(xi) The model will apply to a crash which takes place on a 
sloping road, provided the slope doesn’t vary. Under these 
conditions a skidding car will have constant acceleration, 
and the magnitude of the acceleration will depend only on 
the coefficient of kinetic friction and the slope of the road. 


The force diagram for a downhill skid is shown in Figure 3. 


direction of motion 


Figure 3 
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‘The equations of motion are 

ma = mgsin6 — F, 

N = mgcos6, 

F=y'N, 
where @ is the angle of the slope. These equations lead to 

a = (sin — p' cos)g- 
As before, the acceleration is constant and the same for the 
test car as for the crashed car, provided that ju’ is the same 
for each and that the test is carried out on road with the 
same slope as the road where the crash took place. (We must 
suppose that tan 6 < yi’, otherwise the car will never stop!) 
The model will apply to a Rolls Royce provided that the 
test car has similar tyres. The size and weight of the crashed 
car (let alone the price) are not relevant, so long as the 
coefficient of kinetic friction can be duplicated in the test. 


A crash into a headwind will not be covered by the model, 
because the headwind will introduce an additional force, of 
air resistance, which will depend on the speed of the car 
relative to the wind; so the assumption of constant 
acceleration will no longer apply. 


If the road is wet at the time of the crash, but dries before 
the test, the model will not apply because the assumption 
that the acceleration is the same in the test and the crash 
will not be valid: the coefficient of kinetic friction for a wet 
road will be different from that for a dry road. 


(xii) The formula u = /—2as does predict that u is an 
increasing function of s, as expected. However, by itself it is 
not enough to solve the problem, because we have no direct: 
way of finding the value of a. The test gets round this 
difficulty, and provides a value of a, 80 long as the 
assumptions that have just been discussed are satisfied, 


(xiii) The symbol s is being used with two different 
meanings: as the distance travelled in general by a point 
moving with constant acceleration, and as the length of the 
skid marks of the crashed car. The method used in the 
passage to analyse the skids is to take the general formula 
for motion with constant acceleration in the form 

v® = u? + 2as and substitute the particular values of v, u, a 
and s appropriate to the particular circumstances. In the 
case of the crashed car we must substitute the length of the 
skid marks of the crash for s; since the same symbol! has 
been used for this as for the general distance, we are led to 
the odd looking equation s = s. The moral of this story is 
that it is better to observe the rule that each symbol should 
have only one meaning. 


(xiv) If the skid marks of the crashed car are four times as 
long as those of the test car, then s/s, =4, so u= 2, 
which means that the crashed car was travelling at 80mph 
when it began to skid, and was exceeding the 70 mph speed 
limit. 

(xv) The instructions to a policeman going to the scene of 
an accident for the first time might take the following form. 


1. If there has been a skid which has left a mark on the 
road, measure the length of this mark. 


2. Check that the tyres on the crashed car are in a 
comparable state of wear to those on the patrol car, and 
that the circumstances have not changed since the 
accident in a way which will affect the slipperiness of 
the road (for example, that the weather conditions have 
not changed). 


3. With a clear road, driving the patrol car, approach (at a 
safe speed) the spot where the accident took place, and 
induce a skid. Note the initial speed of the patrol car, 
and measure the length of its skid. 


4. To calculate the speed of the crashed car when it began 
to skid, divide the length of the crash skid by the length 
of the test skid, take the square root of the result, and 
multiply it by the speed of the patrol car (in mph). The 
answer will be an estimate of the speed, in mph, of the 
other car when it started skidding. 


(xvi) The model has been based on the assumption that 
uv =0. If this is not the case, then, provided that the test 
still gives the correct value for the acceleration, we have 


u= Vv? — 2as > V—2as. 


So the actual speed of the crashed car would have been 
greater than that estimated using the model. 


Solutions to the exercises in Section 3 


1. The actual problem described in the summary of the 
model of double glazing is to find how the rate of transfer of 
heat through a double-glazed window depends on the 
difference between the air temperatures on either side of it. 
For many applications it would also be of interest to find 
how the rate of transfer of heat depends on the design of the 
window ~ for example, on the width of the space between the 
two panes, or the thickness of glass used. This information 
is contained in the model, but it is necessary to know how 
the U value depends on such features to obtain it explicitly. 
The most likely use of a model of heat transfer through 
double glazing is in calculations concerned with saving 
energy. For example, to calculate the savings resulting from 
installing double glazing in one’s home, one would have to 
estimate the annual heat loss through the windows with and 
without double glazing; the difference would be the amount 
of energy saved per year, and the cost of providing this 
energy would be the financial saving. This would have to be 
compared with the cost of installing the double glazing. 
Again, any attempt to estimate the reduction in national 
demand for energy resulting from a policy of improving the 
insulation of buildings would have to take into account the 
effects of double glazing, among the various forms of 
insulation. 


2. The shot-putter’s objective is to put the shot as far as 
possible. How far he or she can put it depends in part on 
physique and strength, in part on exactly how the shot is 
thrown. A mathematical model of the motion of the shot 
cannot help much with the former, but it can help with the 
latter. One important factor affecting the range of the shot 
is the angle at which it is launched. The formula relating 
the range to the angle of launch can be used to find which 
angle gives the maximum range. The shot-putter can train 
to use this launch angle, and so try to ensure that he or she 
achieves the maximum range. 


3. The sandwich shop manager will want to order just 
sufficiently many sandwiches to satisfy one day’s demand; 
ordering too many sandwiches will lead to waste, while 
ordering too few will result in disappointed customers. The 
problem therefore is to predict demand. We must next ask 
what factors will affect the demand for sandwiches on any 
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particular day. There are several factors which will affect the 
overall level of demand, such as the location of the shop; 
however, the most important question, so far as the manager 
is concerned, is how the demand varies from day to day. 
There are three obvious factors affecting the daily variation 
in demand: what day of the week it is, what time of the year 
it is, and what the weather is like. So the basic problem is 
to find out how the daily demand for sandwiches will depend 
on the day of the week, the time of year and the weather. 


The distinctive feature of sandwiches is their short shelf-life. 
Another commodity with the same feature is newspapers. 
‘The demand for newspapers will depend on different factors; 
one major factor will be what newsworthy events are taking 
place, if any. 


4. The basic problem underlying the analysis of the time of 
fall of objects dropped from Clifton Suspension Bridge is to 
find how the position of a falling object (modelled as a 
particle) depends on how long a time it has been falling. 
We begin by selecting an origin and an axis, and we use the 
particle's coordinate to define its position. In this case the 
most obvious choice of origin is the point from which the 
particle is dropped, and the most obvious axis points 
vertically downwards, It is also sensible to measure time 
from the moment the particle is released, The main 
variables are therefore 


e x, the particle's coordinate with respect. to the chosen 
axis (x is the distance the particle has fallen, in metres), 


et, the time in seconds since the particle was released. 


It is tempting to consider the height of the bridge as a 
‘variable’. We have used quotation marks here to point out 
the oddity of this suggestion: the effects of any significant 
variation in the height of a bridge do not bear thinking 
about. By widening the problem to cover the whole period 
of the motion, and by using the particle’s coordinate as a 
variable, we avoid this absurdity. 

Position and time are very often the main variables in 
mechanics problems. 


5. In modelling the repayment of a mortgage, the basic 
problem is to find how the amount of money still owed 
depends on how long one has been repaying the mortgage. 
The main variables are therefore 


@ u,, the amount of money owed after r years, in £, 


er, the number of years since the mortgage was taken 
out, 


6. This problem is concerned with deciding the position of a 
supermarket: so one (or more) of the variables must specify 
the position. One way of doing this which would be 
appropriate to this kind of problem would be to use a map 
reference; in effect, a map reference is a pair of coordinates, 
so we could say that the coordinates of the possible sites of 
the supermarket are a pair of main variables. The next, 
question is to ask what determines where the supermarket 
will be sited. The basic factor is presumably the number of 
customers it will attract. Leaving aside such matters as 
whether people prefer one supermarket chain to another, the 
simplest assumption is that people who use supermarkets 
usually shop in the nearest one. For any possible site of the 
new supermarket, there is an area surrounding the site with 
the property that, if you live in this area, and if the new 
supermarket were built at that site, the supermarket closest 
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to where you live would be the new one. This area we call 
the catchment area of the site. The crucial question is how 
many people live in the catchment area. The management of 
the supermarket chain will probably consider only sites for 
new supermarkets for which the catchment area contains a 
certain minimum number of people; and it will seek to 
choose the site to maximize the number of people living in 
the catchment area, The main variables are therefore 


© the number of people living in the catchment area of 
any possible site of the new supermarket, 


¢ the coordinates of each possible site. 


7. The main simplifying assumptions are 


the object may be treated as a particle (as in the 
solution to Exercise 4 of this section), 


¢ the particle is dropped vertically and moves in a vertical 
straight line, 


© (for a first model, at least) air resistance can be ignored, 
so that the only force acting on the particle is its weight, 


8. The fundamental assumption is that the system (the 
double-glazed window and its surroundings) is in a steady 
state: that is, neither the internal or external temperature, 
nor the rate of heat transfer, is changing in time. In addition 
we are assuming that the heat is transferred in one direction 
only (perpendicularly to the plane of the window), so that 
the effects of the edge of the window are ignored; and we are 
ignoring the effects of frames and glazing bars. The window 
is assumed to be closed and draught free. We are also 
assuming standard models of conduction and convection. 


9, For a first model it would be wise to make as many 
simplifications as possible, consistent with retaining some 
degree of realism. Take the windscreen to be flat and 
rectangular. Deal with a single wiper only to begin with 
(that is, solve the problem for the rear window first, or 
assume the car is a Citroén). Assume that the area swept 
out by the wiper blade has the shape shown in Figure 1. 


Figure 1 


‘The curves are sectors of circles with the same centre, and 
the straight sides are segments of radii of the larger circle. 
Such a shape should probably be called a sector of an 
annulus. Finally, assume that the wiper blade never lifts off 
the glass, so that the region swept out lies entirely within 
the rectangular area of the windscreen. 

(The problem we are left with after these simplifications is 
still quite hard enough: it is to find the largest sector of an 
annulus which can be fitted inside a given rectangle, A full 
solution would require, in effect, the specification of the 
dimensions of the wiper: length of arm, length of blade and 
angle of sweep; and the location of the point about which it 
rotates.) 
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10. In each year the building society (or whoever is the 
lender) will charge interest on the amount still owing at the 
beginning of the year; this will be added to the outstanding 
amount. During the year, twelve equal monthly repayments 
are made. When these are subtracted in total from the 
amount outstanding plus interest, what remains is the 
amount owed at the beginning of the next year. 


If the amount owed at the beginning of year r is u, and the 
interest rate is 1% per annum, then the interest. to be added 
is Ju, /100, while the total repayment in the year is 12M if 
the monthly repayment is M. Thus 

Ur = (1 + ) ur — 12M. 
In the example in Unit 1, the interest rate is 15%, so the 
formula becomes 


Ur¢1 = 1L.15ue — 12M. 


(All variables representing sums of money are in units of £, 
whereas J is an annual percentage.) 


11. The force diagram is shown below. 


Figure 2 


The forces acting on the car body are the tension in the 
spring, T, and the resistance due to the damper, R. 
In order to derive the equation of motion, we must first 
define the coordinate. We shall take the x-axis to point 
upwards, with the origin at the equilibrium position of the 
car body (there are other possibilities, of course). The 
equation of motion is obtained by using Newton's second 
law; in terms of the symbols for the forces shown in the 
figure we have 

mi = —R—mg—T. 
We must next express R and T in terms of x and #. The 
resistance force, assumed to be due to linear damping, is 
given by R = ré, where r is a constant. The expression for 
the tension may be derived quickly as follows. For a perfect 
spring the tension will be a linear function of x, In other 
words, we must have T = ke +1, where k is the stiffness of 
the spring and / is an as yet unknown constant. We may 
find | by noting that T+ mg = 0 when x = 0, which 
corresponds to equilibrium. Thus | = —mg, and 
T = kx — mg. The equation of motion is therefore 

mi = —r — mg — (kx — mg), 
which may be rewritten 

mi +ri+kr=0. 


12. To express t in terms of x we must work out the time 
taken to find a parking space, and the time taken to walk 
from there to the shops, and add the two together. The time 
taken to find a parking space is assumed to be inversely 
proportional to x, so it can be written k/x where k is a 
constant. I then have to walk a distance x at a speed of 
(say) v, which takes a time x/v. The total time t is therefore 
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given by 
kw 
t==4=, 
zu 


All quantities here are measured in appropriate SI units — 
or, indeed, in any consistent system of units (yards and 
minutes might be considered more suitable than metres and 
seconds). The constant of proportionality & must have units 
of ms (in the SI system). 

‘The question asks only for the formulation of the problem, 
but it is easy to find the value of x which minimizes t, and it 
turns out to be Vkv. (Note that the units make sense: since 
k is in ms and v in ms~', the result is in metres as 
required.) I know my own walking speed, or can easily find 
it; the difficulty with this model will be to estimate k, One 
reason for this is that although it is stated that k is a 
constant, it will actually vary from day to day, and even 
during the day. In fact, the value of k will depend on how 
much competition there is for parking places, that is, on how 
many other people want to go shopping at the same time as 
Ido, 


13. The solution to the logistic equation, 

P M 

T+ (M/Po — le" 

has the following features that suggest that it would be a 
satisfactory model of the growth of a population. In the first 
place, P is an increasing function of t (provided that 
Py < M), so it does indeed model population growth. 
(In fact, the same general formula would model a growing, 
stable or declining population, according as Po <M, 
Py = M or Py > M.) Secondly, there is an upper limit to P, 
namely M, and P approaches M more and more closely as 
t— oo. This accords with the expectation that a population 
cannot increase indefinitely, but that there is a natural 
upper limit to its growth, 
Thirdly, if Po (the initial size of the population) is small 
compared with M, then for small values of t, P ~ Pye®*. 
(This result is shown in detail on page 21 of Unit 9.) This 
means that a population which is small initially will grow 
exponentially to begin with. This is the expected growth 
pattern of a population while there is no competition for 
food and space. These features of the solution of the logistic 
equation are immediately apparent from its graph (shown 
below); indeed, sketching the graph would in many ways 
have been the best way of answering the question. 


Population, P 
M. 


Py 


Figure 3 


There are two ways of rearranging the formula to provide a 
way of checking the model against data. The first is to 
rewrite it as 


(-)e"= (3-9), 
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and then take logs to obtain 


log. (#-1) —at = log, (4-1). 


This says that a plot of log,(M/P — 1) against t should give 
a straight line. In order to use this method, however, one 
must first estimate M, the size of the equilibrium 
population; this will not be easy if the population is still 
growing. There is an alternative method, which does not 
have this drawback, which can be used if the data are for 
the population sizes at equal time intervals. Let P,, be the 
population when t = n7', where T is some fixed length of 
time (one year, for example). Then 


(6 = :) = (#-1) canter 


= (# = 1) ete cx (x - 1) et 
from which we obtain . 


Pri Pa _ (l=e™ oT 
Fos Pa - ( a ) Poss = (6 Ei 
It follows that a plot of (Pas; — Pn)/Ps against Posi 
should produce a straight line. (See Subsection 2.4 of Unit 2 
for more details of these methods.) 


(The logistic equation is not a satisfactory model for the 
growth of every population, of course. It takes no account of 
population movements, of the interaction of the population 
with its environment (overgrazing, the effects of predators), 
or of random influences like the weather, to mention just 
three of many factors which could affect population growth, 
Populations of species for whose development these are 
important factors will not follow the logistic pattern of 
growth. Nevertheless, the logistic equation is a very useful 
basic model in population dynamics, much as motion in one 
dimension under gravity is useful in particle dynamics 
despite the factors — like air resistance — it leaves out of 
account.) 


14. The general solution of the recurrence relation for the 
mortgage problem (with interest at 15%) is 

tia = B(1.15)" + 80M. 
‘The initial condition, that uo = 20000, can be used to 
determine the constant B in terms of M, giving 

B = 20000 — 80M. 


We wish to determine the monthly repayment required to 
pay off the loan after 25 years: this leads to the equation 

uzs = 0 = (20000 ~ 80.M)(1.15)*" + 80M, 
or 

20000(1.15)"> 

M= 30(@.15) —1) j = 257.8324... 
(These results are derived in detail on page 36 of Unit 1.) 
The practical problems of implementing the solution arise 
because the calculation of was is ill-conditioned with respect 
to small changes in the value of M; the monthly repayment 
is a sum of money, and so the exact value of M cannot be 
used, but must be rounded. As noted in the description of 
the model, rounding M up to 258, and working throughout 
to the nearest whole number of pounds, results in an 
overpayment of £427. If, on the other hand, one works 
throughout to the nearest penny, taking the monthly 
repayment to be £257.83, then after 25 years the sum of 
£5.99 will still be owing. This is despite the fact that the 
‘error’ in each repayment is 0.24p, which over 25 years 
accumulates to only 72p. 


MST204 M Solutions 


In fact, mortgage calculations are very sensitive to changes 
in almost any of the relevant data. One good illustration of 
this is the fact that when the interest rate is 15% and the 
loan is £20000, a monthly repayment of £250 is required 
simply to pay the interest. In other words, if the monthly 
repayment is £250 then the amount owed stays constant at 
£20000, But if the monthly repayment is £258 the loan is 
more than paid off after 25 years; while if the monthly 
repayment is £249 then the amount owed will actually 
increase. 

The moral of this exercise is that it is often not enough just 
to calculate a numerical answer to a modelling problem. 
One should also consider the sensitivity of the answer to the 
accuracy of the data, and the effects of any approximations 
used. 


15. First, one can check that the formula makes sense. Note 
that t is an increasing function of m, and that t = 0 when 
m=0. The graph of km?/* is shown in Figure 4. 


Figure 4 


‘The formula therefore predicts that as the mass of the piece 
of food to be cooked (a potato, for example) increases, the 
cooking time also increases; but that twice the mass of food 
doesn’t need twice the cooking time, but only about 60% 
longer (27/9 = 3/4 ~ 1.59). These predictions seem 
qualitatively correct. 


< 


‘There are, in principle, two ways of checking the formula in 
more detail. One is to carry out some experiments. The 
other is to find a theoretical derivation of the formula. 


To check the formula by experiment, you would need to 
cook a number of samples of the same food but of different 
masses. In order to tell when the food was cooked you would 
need a cooking thermometer or temperature probe of some 
kind, which could be stuck into the sample being cooked and 
which would register the temperature at its centre. Each 
sample would have to be cooked by itself, the oven 
temperature being the same for each. 


You would record how long it took for the temperature of 
each sample to reach some predetermined value. When the 
experiment was complete, you would plot the logarithm of 
the cooking time against the logarithm of the mass for all 
the samples, and expect to obtain a straight line of slope 3. 


The concepts of Unit 12 could be used to derive the formula 
theoretically, though the problem is much more complicated 
than any considered in that unit. Even so, we can note that 
the main heat transfer process involved in the heating of the 
food is conduction, which suggests that the value of the 
constant of proportionality k will depend on the thermal 
conductivity of the food, and also the oven temperature. It 
is also worth remarking that the rate of transfer of heat into 
the food will be proportional to the area of its surface; the 
surface area of a sphere is proportional to the square of its 
radius, and therefore to the 2/3 power of its volume, or 


aT 


equally the 2/3 power of its mass (supposing the density to 
be uniform). This is supportive evidence for the formula, 
though not of course a conclusive derivation. However, this 
approach does suggest that the shape of the object being 
cooked is not a significant factor, because the relation 
between surface area and volume just derived holds true, in 
a sense, for any shape. This remark is based on the 
observation that a large baking potato (say) is much the 
same shape as a small potato; the larger one is simply a 
uniformly scaled up version of the smaller one. But when 
you uniformly scale an object by a factor A, you scale its 
surface area by a factor \*, and its volume by a factor \". 
This means that, for any collection of similar objects 
(potatoes, for example), the surface area of any one of them 
is proportional to the 2/3 power of its volume, or the 2/3 
power of its mass. If the hypothesis that cooking time is 
proportional to surface area is correct, then the cooking time 
for one of these objects is indeed proportional to m?/" where 
m is its mass. 

To find how long it would take to bake a potato, one would 
need to know the value of the relevant constant &, This will 
take the same value for all potatoes cooked at a given oven 
temperature, but will depend on the oven temperature 
chosen. The choice of oven temperature has a lot to do with 
the eating quality of the food when it is cooked. You must 
select an oven temperature that produces baked potatoes as 
you like them; it will then be sufficient to carry out one test 
baking to determine the appropriate value of k. Armed with 
this information — which, presumably, you will have recorded 
somewhere for future use — you need only weigh your potato 
to work out the cooking time. In order to avoid 
disappointment, it would be wise when recording the value 
of k to make a note of the system of units you have adopted, 
so that you will know whether to weigh the potato in grams 
or ounces. This is a situation in which SI units are not the 
most useful, since cooking times are rarely measured in 
seconds. 


(Cookery books — that is, books of recipes — tend to avoid 
such sophisticated matters as 2/3 powers. I know of only 
one cookery book in which the 2/3 power formula is 
discussed, and that is an old one: The Grammar of Cookery 
by Philip Harben, Penguin 1965. Harben claims that the 
formula ‘holds good for every time calculation in cookery, be 
it roasting, boiling eggs, baking cakes, deep-frying — every 
single thing’. But he doesn’t expect the reader to use the 
formula to work out the cooking times: he gives a table of 
cooking times in each recipe.) 


16. The exponential model is fairly accurate for many 
populations in a state of rapid increase, over relatively short 
periods of time. Its major drawback is that it predicts 
unrestricted and indefinite growth, which cannot continue in 
fact since resources are finite. 


In the exponential model it is assumed that the rate of 
change in the population is a constant multiple of the 
population size; in other words, that the proportionate 
growth rate of the population is constant. (The 
proportionate growth rate of the population is the difference 
between the proportionate birth rate and the proportionate 
death rate. The assumption that each of these is separately 
constant will lead to the exponential model, but in fact it 
holds in more general circumstances.) 
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To model competition for scarce resources we assume that 
the larger the population is, the smaller its rate of 
reproduction will be, or the larger its mortality rate, or 
both. This amounts to saying that the proportionate growth 
rate of the population is not constant, but depends on the 
population size; and that it is a decreasing function of 
population size. In the absence of any detailed information 
about exactly how the proportionate growth rate g(P) varies 
with the population size P, we make the simplest possible 
assumption, which is that it is a decreasing linear function 
of P: we may thus write 


aP)=a(1-F), 


where a and M are positive constants. This particular form 
is chosen because the constants a and M are relatively easy 
to interpret in terms of the model. Tt follows from the 
assumptions that the population size P must satisfy the 
differential equation 

1dP _ (1 F\ 

Pa °\’" mM)’ 
which can be written 

dP P 

=P (1-3): 
‘This is the logistic equation. We see that M is the value of 
P for which dP/dt = 0, which means that a population of 
this size will remain at the same size indefinitely; so M is the 
equilibrium population size. When P is small the differential 
equation approximates to the exponential equation 


GaP 


this shows that the logistic model incorporates the 
exponential model as an approximation for small 
populations, and that a is the proportionate growth rate 
when the population is going through this approximately 
exponential phase of growth. 


17. In the Clifton Suspension Bridge problem, the result of 
introducing the resistive force, assumed proportional to 
speed, is to change the equation of motion from 


m& =m, 
ia 
to 
dv 
eo mig — eh 
moe = mg — ke. 


where k is a positive constant. The qualitative effect of this 
change on the predicted motion is quite significant. For 
example, the simple model predicts that the speed of the 
falling body will increase linearly with time, and that in 
principle the speed could become arbitrarily large. But 
when the resistive force is included, the speed can never 
exceed the terminal speed vp = mg/k. On solving the 
equation of motion for v in this case, we find that 


y= 7s (1— eM) = up (1 eM). 


(See page 42 of Unit 4 for the full derivation of this result.) 
It follows that when ¢ is large, so that e~*'/" is small, 

uv ~ wp. When t is small, on the other hand, 

e~*t/m ~ 1 — kt/m (using the first two terms in the Taylor 
series of the exponential, which is given in Subsection 6.4 of 
the Handbook); thus 


ve 4 (1 — (1 —kt/m)) = at, 


which agrees with the result in the absence of the resistive 
force. 
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Figure 5 


The predicted time of fall will be longer with the resistive 
force than without, since the speed is smaller. The relation 
between distance fallen and time is obtained by solving the 
differential equation 


de = “ket /m) 
he )i 
the solution turns out to be 
mop ntfs 
oe ee 
Again, it is interesting to consider what happens for large 
and for small t. For large t the term e~**/" is small, and the 


motion is given approximately by 
m 
exer (t-F), 


which is that of a body falling with constant speed vp, For 
small t we have e~"*/" ~ 1 — kt/m + 3(kt/m)?, whence 


2 & upt— Sa (1 = (1 = kt/m + 3(kt/m)*)) 


x=vpt- 


= igt. 
Once again, we obtain the result for the motion under 
gravity alone as an approximation for small ¢ (though we 
have to include an extra term in the Taylor expansion of the 
exponential). 
To find how long it takes to fall a specific distance h, we 
would have to solve the equation 
mi =Kt/m 

h= opt — 2 (1—e™/) 
for t. Unfortunately, it is not possible to solve this equation 
explicitly to give t in terms of h. However, a sketch of the 
graph of x against ¢ shows that the equation has a solution 
lying between /2gh and h/vp +m/k. 


Figure 6 


It is also worth noting that one consequence of introducing 
the resistive force is that the results depend on the particle's 
mass. For motion under gravity alone, this is not so: all 
bodies fall equally fast (as Galileo is supposed to have 
demonstrated in Pisa). But when air resistance takes effect, 
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two bodies of the same dimensions, one made of lead and 
one of gossamer, will have different times of fall. 


18. The effect of taking the drug orally is that the drug in 
the bloodstream will be continually replenished as what 
remains in the stomach is absorbed. This is similar to 
intravenous infusion, but with one important difference: the 
stomach will not pass the drug to the bloodstream at a 
constant rate. It seems likely that the absorption of the drug 
in the stomach works something like the elimination process 
in the original model: the rate at which the drug is absorbed 
(and therefore transferred to the bloodstream) is 
proportional to the amount remaining in the stomach. This 
implies that the amount of drug in the stomach decreases 
exponentially, so we can write 
Q(t) = De, 
where Q(t) is the quantity of drug in the stomach, in mg, t 
seconds after the dose was taken. Here D = Q(0) is the 
initial dose, again in mg, and y is a positive constant which 
measures the rate of absorption. Now drug lost from the 
stomach is gained by the bloodstream, so that the rate of 
gain of drug by the bloodstream is —dQ/dt = »De~“'. The 
equation for the quantity of drug in the bloodstream 
becomes 
a =-Aq+pDe~""', 
the first term on the right-hand side representing excretion 
via the kidneys (as before). 
This equation can be solved readily by the integrating factor 
method (for example). The general solution, for \ # 4, is 
= Aew™t 4 HD o-ut, 
OT i earn 
the solution which satisfies the initial condition q(0) = 0 is 
uD 
A~ 
‘The drug concentration is therefore 


re (et eM), 


mat gay, 


where V is the apparent volume of distribution. 
‘The graph of this function looks like Figure 7. 


Figure 7 


‘The concentration is initially zero, grows to begin with at 
the rate 4D/V, reaches a maximum when 


p= 108A = log. w 
eas 


and then tails off towards zero again. 

Before we could make any quantitative statements about the 
size and frequency of the dose in this case, we would have to 
know the value of the constant yz, which is a measure of the 
rate at which the drug is transferred from the stomach to 
the bloodstream. However, we can see how to proceed once 
this is known. The initial dose D must be chosen so that the 
maximum value of the drug concentration lies below the 
level at which toxicity occurs. It will take some time before 
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the concentration is sufficient for there to be a therapeutic 
effect; and second and subsequent doses must be 
administered before the concentration falls below the 
therapeutic range, in order to allow for this delay. 


Solutions to the exercises in Section 4 


1. The process to be modelled is the warming of something 
cold — a liquid contained in a bottle - when it is placed in 
surroundings which are at a higher temperature than it is. 
The problem, therefore, is to find how the temperature of 
the liquid varies with time. Let ¢ be the time, in seconds, 
since the bottle was placed in the warmer surroundings; and 
let A(t) be the temperature of the liquid, in degrees 
Centigrade, at time t. One would expect @ to be an 
increasing function of t, and to approach the temperature of 
the surroundings more and more closely as t increases but 
never to quite reach it. 

The liquid’s temperature increases because heat energy is 
continually being transferred to it from its surroundings; the 
most significant mechanism of heat transfer in this case is 
likely to be conduction through the wall of the bottle. To 
model this it will be necessary to make some assumption 
about the shape of the bottle: the simplest assumption 
which is at all realistic (in most cases) is that the bottle can 
be regarded as a cylinder. The heat conducted through the 
wall of the bottle will in fact warm both the liquid and the 
bottle; to simplify matters, the warming of the bottle itself 
will be ignored, so it will be assumed that all of the heat 
warms the liquid. It will also be assumed that the 
temperature of the surroundings remains constant as the 
liquid warms up. The effects of convection will be ignored, 
as will the opening of the top of the bottle. 

In the steady state, the rate of transfer of heat by conduction 
through the side and ends of a cylinder of length I metres, 
radius r metres and thickness b metres is q watts, where 


2nl mr? 
= ——— + 2— } AO, 
, { iog.((r+8)/) * } 
where x is the thermal conductivity, in Wm~! °C7!, of the 
material of which the bottle is made (glass or plastic), and 
46°C is the difference in temperature between the inside 
and the outside of the bottle. Put 
2nl 
ot SE op 
"\log.(r+8)/r) +" 6 } ze 
to save writing. Then in a short time 6t the liquid gains an 
amount of heat energy 
git = K(0, — @)ét, 
where @,°C is the temperature of the surroundings. The 
resulting increase in the temperature of the liquid, 60°C, is 
given approximately by 
RN Cie 
me me 
where m kg is the mass of the liquid and cJ kg~' °C? its 
specific heat. On letting 6t tend to zero we obtain 


— =—(6,— 6). 
me 
This differential equation is the required model of the 
warming process. 
If we put ¢ = @—@, and A = K/(me), the equation becomes 
do 


= — >, 
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as in the text. The solution is 

= H%e™, 
or 

O(t) = 0, + (80 — O.)e™*, 
where 4 is the initial temperature of the liquid. Note that 
since 0) < @,, the temperature of the liquid is in fact 
increasing and tending to 8, from below: perhaps it is better 
to write the solution as 

4(t) = 0, — (0. —O%)e™, 
to make this a little more obvious. 
Suppose that the temperature at which the liquid becomes 
unpleasant to drink is 6, °C. The time it takes for the liquid 
to reach this temperature is then 

1 1 (2 = 9% ) 

X NG, Oy 
To find a numerical value, the values of the various 
parameters need to be known. 


The results of this analysis are very similar in appearance to 
those of the model of the cooling of tea considered in the 
text, even though that is a problem of cooling by convection 
while this is a problem of warming by conduction. In fact, 
any situation involving heat transfer processes which, in the 
steady state, lead to a rate of heat transfer which is 
proportional to temperature difference, will give a similar 
answer. Of course, the nature of the heat transfer processes, 
the shape of the container and so on will affect the outcome 
~ but they will do so only through the value of \. This 
observation suggests a different approach to the evaluation 
of parameters. One could ignore the details of the process, 
and find the value of \ directly by experiment instead. Thus 
to deal with the warming of a bottle of milk, it would be 
enough to conduct one experiment to find the value of \ for 
a milk bottle; this could be done by noting that the time 
taken for the temperature difference between the contents of 
the bottle and the surroundings to halve is (log, 2)/A. With 
the value of A known, the formula derived above will give 
the time it takes for a bottle of milk to reach any given 
temperature when taken from the refrigerator, and it can 
also be used to determine what the temperature of the milk 
will be a given time after it has been collected from the 
doorstep and put into the refrigerator. 


2. ‘The temperature inside the house while it is cooling 
(with the heating off) obeys the law obtained in the text 
and the previous exercise, namely 

O(t) = Ox + (Bo — BxJe™, 
where now 4(t) °C is the temperature in the house t seconds 
after the thermostat has turned the heating off. @, °C is the 
temperature outside the house (assumed constant), and 
@o °C is the temperature at which the heating cuts out 
(called the ‘predetermined higher value’ in the question). 


Here ) is a parameter which measures the rate of heat loss 
for the whole house, which could be defined as (log, 2)/T 
where T is the time taken for the temperature difference 
between the inside and the outside of the house to fall to half 
its initial value: see the solution to the previous exercise. 
‘The heating remains off until the temperature inside the 
house falls to the ‘predetermined lower value’, at which 
point the heating comes on. By assumption, the heating 
system supplies heat at a constant rate, say H watts. But 
while the heating is on, the house continues to lose heat by 
the same processes as before, as well as gaining heat from 
the heating system, We must therefore modify the model of 
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cooling to take account of the heat supplied. Formerly, we 
had the expression 

gét = K(0, — 0)ét, 
for the amount of heat transferred in time ét; here 
K = \me, where now mkg is the mass and cJ kg~! °C~ the 
specific heat of the air in the house. To take account of the 
heat supplied we must change this to 

qét = {K(0, — 0) + Hot. 
The differential equation for the temperature inside the 
house now becomes 

dO H 

a = NO - 8) + ; 
Before rushing in to solve this equation as it stands, let us 
note that it predicts that d@/dt = 0 when 

H 

0=6,+ rerd 
‘Thus the temperature of the house can stay constant at this 
value, with the heating on. We call this temperature the 
equilibrium temperature and denote it by @,. At the 
equilibrium temperature the heat supplied just balances the 
heat lost to the outside. In terms of @, the differential 
equation can be rewritten as 


7) 
T= NGn ~ 8). 


‘The solution has the same form as before, but with 8, 
replacing @,; so in the heating phase we have 

O(t) =O. + (0:1 — OnJe™™, 
where now @(t) °C is the temperature in the house t seconds 
after the thermostat has turned the heating on, @, °C is the 
equilibrium temperature when the heating is on, and 0 °C is 
the temperature at which the heating cuts in (called the 
“predetermined lower value’ in the question). We must 
assume that the heating system is sufficiently powerful that 
6. > 4, since otherwise the temperature inside the house 
will never reach the value at which the heating cuts out, 
Assuming this to be so, the heating will continue to run 
until the temperature 4 is reached, when the thermostat 
will operate again, and the cooling half of the cycle will 
recommence, The temperature in the house will vary as 
shown in Figure 1. 


Figure 1 


The difference between the equilibrium temperature 0, and 
the external temperature 4, depends just on the rate of 
supply of heat and the thermal characteristics of the house; 
so when the external temperature drops, the equilibrium 
temperature drops by the same amount. But the thermostat 
temperatures 4 and 6; do not change. So when it is colder 
outside, the cooling and heating curves move downwards 
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relative to 0 and 6;, with the result that the cooling curve is 
steeper between 4 and @;, and the heating curve less steep. 
Thus when it is colder outside, the heating part of the cycle 
is longer and the cooling part shorter, as one would expect. 
This problem has some features in common with the drug 
therapy model: the range of temperatures allowed by the 
thermostat is like the therapeutic range for the 
concentration of the drug; while the heating part of the 
cycle and the intravenous infusion method of administering 
a drug are modelled in very similar ways. Models occurring 
in apparently totally different circumstances can be 
surprisingly similar; it can often be helpful when modelling 
something unfamiliar to ask yourself ‘have I ever come 
across anything similar before?’. 


Solution to the exercise in Section 5 


1. (i) Things like convenience of handling and stacking 
must, presumably, be taken into account; but surely the 
most important thing, from the manufacturer's point of 
view, is to use as little tin-plate as possible consistent with 
the fact that the can has to hold a given amount of beans 
(or whatever foodstuff it is destined for), 


(ii) The problem is: to find a formula for the area of 
tin-plate required to make a cylindrical can which will 
contain a specified volume of food; this formula should give 
the area in terms of the dimensions of the can (its height 
and the radius of either end); and then to use this formula 
to find the dimensions of the can for which the area of 
tin-plate used is a minimum, 

Tt seems that canned food is sold by volume rather than 
weight (or more properly mass). The evidence for this comes 
from my kitchen cupboard, where I find cans of different 
foods the masses of whose contents range from 397g to 
432g; the cans are identical in size and shape. It is clearly 
better to have a standard size and shape of can, and 
therefore a standard volume, rather then slightly different 
sizes of can for 400g of peas, 400g of beans and 400g of 
pears (for example). 


(iii) The fundamental variables are the area A of tin-plate 
used in the manufacture of a can, the height A of the can, 
and the radius of an end, r, In SI units the linear 
dimensions would be in metres; but in this case the metre is 
rather too large a unit, the millimetre rather too small, and 
centimetres would seem to be the better choice. Then A 
must be measured in cm. Of course, one could use inches 
and square inches. 


(iv) The volume is a parameter rather than a variable. The 
fact that we are concerned only with cans which enclose a 
given volume will enable us to formulate a relation which 
must be satisfied by the dimensions of the can, which will 
enable us to eliminate one of them (either the height of the 
can or the radius of an end) and express the area entirely in 
terms of the other. 


(v) The can is made up of one rectangle of sides h and 2nr 
(which forms the curved body of the can) and two circular 
ends, each of area xr*, Thus 


A=2nrh + 2nr?. 

The volume V which the can encloses is 
V=mrh. 

Eliminating h gives 


A=2nr (4) 42m? =2 Y sar). 
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(vi) The area of tin-plate required to make the seams has 
been ignored. So also has the waste resulting from cutting 
circular pieces for the ends of the can from sheet material. 


(vii) We have A — oo as r + 0 and as r — oo. In the 
former case the can becomes long and thin; in the latter, 
short and fat. 


(viii) Since the area is always positive, it must have at. least 
one local minimum. In fact, it is easy to sketch the graph of 
A by noting that it is the sum of a quadratic function and a 
multiple of 1/r: the graph looks like this. 


Figure 1 


There would appear to be just one local minimum, which 
will provide the solution to the problem. 


We have 


dr a 
and so dA/dr = 0 where 
Vv 
22 
wae 


It is clear from the graph that this is a minimum; but to 
confirm, note that 
@A Vv 
o*=4 ( + ) i 
which is positive for all positive r. 
The ratio of height to radius of a can of minimum area is 
given by 

Vv 
teks 
This means that for a can of minimum area, the height is 
equal to the diameter of the base; in other words, seen from 
the side the can will appear to be square. In fact, most cans 
are taller than they are wide. 


(ix) The most significant source of error is the fact that we 
have ignored the wastage in cutting the circular ends; in 
computing the area of tin-plate used to make the can, we 
should include the waste as well as the tin-plate which 
actually ends up as the can. 


(x) The simplest assumption is that each end is cut from a 
square of material of side 2r, so that the formula for the 
area becomes 


he at 
r 
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The minimum now occurs where 


—2¥ 4167 =0, 
. 


¥ 


8 
‘The ratio of height to radius becomes 
v 
Ue aE os 
7 or 
This gives a shape more like that of actual cans, but still 
rather too short and fat: the corresponding ratio for an 
actual can is 3.10. We still have to take the seams into 
account, but it is hard to believe that this will make a 
substantial difference. This analysis suggests that there is 
money to be saved in changing the shape of tin cans. 
(xi) For a square section can of height h and width w the 
area of tin-plate required is 
A=4hw + 2w? 
(there is no waste). The volume is 
V =hw*, 
so we may write 
aa eae, 
w 
The minimum occurs where 
AY tum, 
w 
or 
w=V9=h; 
the optimal shape is a cube (as one might have guessed), 
and the area required is 6V?/*, For the solution found in the 
previous part we have r= 3V/*, so 
A=4v?/5 +2V?/9 = 6V2/9 also. Thus square section cans 
use the same area of tin-plate for a given volume. 
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